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DESCRIPTION 

DEVICE FOR MEASURING DYNAMIC MATRIX SENSITIVITY 
OF INERTIA SENSOR AND METHOD FOR THE MEASUREMENT 

Technical Field: 

This invention relates to a device for measuring the dynamic matrix sensitivity of an 
inertia sensor and a method for the measurement. Particularly the invention relates to a 
device for the measurement of the dynamic matrix sensitivity of an inertia sensor for a 
varying use such as, for example, an inertia sensor relating to inertial navigation devices to 
be mounted on automobiles, submarines, missiles, and airplanes, an inertia sensor to be 
used for motion control of robots, an inertia sensor to be used for measuring the motion of 
a human body, the vibration exerted on a human body, and the motion of an animal, and an 
inertia sensor to be used for preventing image devices and screen images from blurring and 
a method for effecting the measurement. 

Background Art: 

As one example of the inertia sensor, the accelerometer which is well known may 
be cited. It is generally a one-axis accelerometer which is furnished with one sensitive axis. 
When this one-axis accelerometer is calibrated, the calibration is effected by causing the 
direction of motion generated by a motion generating machine to coincide with the axis of 
sensitivity. The degree of freedom of motion to be used for calibrating accelerometers, 
therefore, is a single degree of freedom. The primary calibration using a laser 
interferometer which is reputed to have the highest precision also uses this technique. 

However, since the device which calibrates the one-axis accelerometer mentioned 
above generally entails a three-dimensional motion, it is rare that the device will be limited 
to a one-dimensional motion. The fact that the calibration is carried out by causing the 
direction of motion generated by the motion generating machine to coincide with the axis 
of sensitivity as described above occurs in the calibration which resorts to the measurement 
of the amplitude of the acceleration when the direction of motion is known in advance. 




As concrete examples of the conventional one-axis accelerometer, a piezoelectric 
type accelerometer, an electromagnetic type servo accelerometer, an interference type 
optical fiber accelerometer and a strain gauge type accelerometer have been known. 
Owing to their structures and the natures of their materials, accelerometers are influenced 
by the acceleration components not parallel to the sensitivity axis when the directions of 
application of acceleration to the acceleration sensors fail to coincide with the direction of 
sensitivity axis. 

It is, therefore, apparent that concerning practical motions, the calibration technique 
alone in the present state of affairs has not fully satisfactorily established a method for 
evaluating the performance of an acceleration sensor or perfected a measurement standard 
for the determination of acceleration. 

It is derived that in terms of vector space with three-dimensional transverse motion, 
even the cross or transverse sensitivities of one-axis accelerometers are expressed by two 
parameters as will be subsequently explained. The practice of denoting the two kinds of 
transverse sensitivity by and Sz, y and designating as a magnitude of not more than 
5% and Sz, y as a magnitude of not more than 3% has never been in vogue heretofore. 

It is natural that the one-axis accelerometer generally emits an output signal in 
response to an input component in the direction of the sensitivity axis thereof It is also 
characterized by emitting output signals in response to input acceleration components from 
two directions perpendicular to the sensitivity axis thereof The reason for this property is 
that the piezoelectric accelerometer, the electromagnetic servo accelerometer, the 
interferometer type optical fiber accelerometer, or the strain gauge type accelerometer 
mentioned above is provided with a mass capable of also moving, though slightly, in a 
direction other than the direction of sensitivity axis or something equivalent thereto and, 
therefore, is so configured as to detect the relative motion of this mass or detect a voltage 
or an electric current necessary for preventing this relative motion. 

Heretofore, the accelerometer is set on the one-axis motion generating machine and 
the sensitivity axis of the accelerometer is caused to coincide with the direction of the 
motion generated by a motion generating machine as illustrated in Fig. 3. The concept of 
enabling the accelerometer to be calibrated most accurately by measuring a motion with a 
laser interferometer under such set conditions as mentioned above and consequently 
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establishing a standard for measurement of acceleration is officially approved by the Treaty 
of the Meter as well. Generally, the reference accelerometer is calibrated in accordance 
with the method embodying this concept. 

Then, in industries, it is supposed to calibrate a given accelerometer based on a 
reference accelerometer mentioned above by joining in series connection the reference 
accelerometer which has undergone measurement by the method of Fig. 3 and the given 
accelerometer as illustrated in Fig. 4(b), causing the sensitivity axis to coincide with the 
direction of motion generated by a motion generating machine, and comparing the output 
signals from the two accelerometers. 

The conventional method of calibration which resides in determining the transverse 
sensitivity from the output signal due to a motion only in one direction perpendicular to the 
sensitivity axis as illustrated in Fig. 4(a) and Fig. 4(b), however, is essentially in error in the 
elementary sense. In the sense that this method is an expedient and is capable of 
determining only one transverse- sensitivity, the thought directed toward the decomposition 
and the synthesis of vector supports a judgment that this method views the phenomenon 
only in a two-dimensional space. 

The transverse sensitivity is determined by imparting vibration only in one direction 
perpendicular to the sensitivity axis as illustrated in Fig. 4(c). 

For the sake of surveying the transverse sensitivity more specifically, the behavior 
of the piezoelectric type accelerometer using a piezoelectric material, for example, will be 
explained below. The piezoelectric type accelerometer possesses transverse sensitivity 
because the piezoelectric constant comprises a shear component. That is, the 
piezoelectric substance generates an electric charge which transmits a signal via an 
electrode even to slippage. Generally, in the region in which the voltage (or electric 
current) generated in response to an input signal (acceleration) possesses linearity, the 
sensor sensitivity is defined by the ratio of their magnitudes. Thus, the following formula 
is established. 

The sensitivity axis output voltage (ao X (co)exp(jcot)) of accelerometer 

= normal sensitivity x input component of acceleration in normal sensitivity 
direction + cross (transverse) sensitivity 1 x input component of 
acceleration in direction 1 perpendicular to normal + cross 



(transverse) sensitivity 2 x input component of acceleration in 
direction 2 perpendicular to main axis 
= Sx^(co)aixexp(jcot) + Sx, y ai y exp(jcDt) + S^a^p^o*) 
When this formula is rewritten in the matrix form, the formula 1 is obtained. Here, 
the amplitude of the vector acceleration exerted on the accelerometer is denoted by (a*, a iy , 
aiz) and the time change component by exp(jcot). 
(Mathematical 1) 
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The drawing of the acceleration vector A applied to the accelerometer and the 
decomposition of the vector in X, Y and Z axes is shown in Fig. 2. 

It is well known that the acceleration is a vector which is expressed by amplitude 
and direction. Further, in order that the accelerometer may correctly measure the 
acceleration, the accelerometer must be calibrated with the acceleration as a vector. The 
conventional method of calibration, however, effects the calibration with the magnitude of 
amplitude because the direction of the exerted acceleration is determined at the stage of 
setup. 

When the accelerometer is put to actual services, there are times when the direction 
of motion can be forecast and there are times when the forecast cannot be attained because 
of the possibility of the accelerator being exposed to application of acceleration in every 
direction. 

In the case of an earthquake or accidental collision of cars, it is not possible to 
know in advance the direction of motion. When the calibration is made only in amplitude 
(one-dimensionally) as practiced conventionally, there are times when no correct magnitude 
of acceleration can be obtained. Thus, the desirability of calibrating the accelerometer by 
applying acceleration to the accelerometer in all actually conceivable directions has been 
finding general recognition. 
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The technical background of this invention has been described heretofore. As 
concrete examples of the prior art of this invention, the following documents have been 
known. 

[Non-Patent Document 1] Vibration Engineering Handbook, complied by Osamu 
Taniguchi, published in 1976 by Youkendo, Chapter 13 "Determination of Vibration," 
13.3.2 "Calibration of vibration measuring device" 

[Non-Patent Document 2] ISO (the International Organization for Standardization) 
16063-11: 1999(E) 

[Non-Patent Document 3] FINAL REPORT ON KEY COMPARISON CCAUV. 
V-Kl Hans-Jurgen von Martens, Clemens Elster, Alfred Link, Angelika Taubner, 
Wolfgang Wabinski, PTB-1.22 Braunschweig, October 1, 2002 

[Non-Patent Document 4] ISO 5347 part 11 Testing of transverse vibration 
sensitivity 

[Non-Patent Document 5] ISO 5347 part 12 Testing of transverse shock sensitivity 
[Non-Patent Document 6] ISO 8041 Human response to vibration - Measuring 
instrumentation 

[Non-Patent Document 7] ISO 2631-1, 1997 Evaluation of human exposure to 
whole-body vibration Part 1 : General requirement 

[Non-Patent Document 8] ISO 5349-1, 2001 Measurement and evaluation of 
human exposure to hand-transmitted vibration - Part 1 : General guidelines 

Disclosure of the Invention: 

This invention relates to a device for measuring dynamic matrix sensitivity of an 
inertia sensor and serving to enable an inertia sensor to be calibrated by application thereto 
of a vibration possessing an acceleration component along the degree of not less than two 
degrees of freedom selected from among the degrees of six degrees of freedom of motion 
and a method for the measurement. The first aspect of this invention is characterized by 
being provided with a motion generating machine capable of inducing a translational 
motion or a rotary motion, an acceleration measuring device which as a device subject to 
calibration is at least temporarily fixed on the motion generating machine, an angular 
velocity measuring device or an angular acceleration measuring device, an output means for 
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fetching an output from the device subject to calibration, a single or a plurality of light 
reflecting materials, a displacement measuring means capable of seizing a multidimensional 
motion by the use of a laser interferometer formed by radiating laser beams from a plurality 
of directions to the single or plurality of light reflecting materials, a data processing device 
for processing the data indicating a multidimensional state of motion obtained from the 
displacement measuring means and converting the resultant multidimensional translational 
motion or multidimensional rotary motion into magnitudes fit for a predetermined 
coordinate system, and a display means for displaying or a conveying means for conveying 
the output of the data processing device and the output of an acceleration measuring device, 
an angular velocity measuring device, or an angular acceleration measuring device which is 
the device subject to calibration. 

The second aspect of this invention is characterized by the fact that the motion 
generating machine mentioned above generates a periodic motion. The term "periodic" as 
used herein means that the vibration produced in each cycle by the motion generating 
machine mentioned above is periodic within the period in which no influence is exerted on 
the measurement using the vibration of the next cycle. 

The third aspect of the invention is characterized by the fact that the motion 
generating machine mentioned above generates a motion of the nature of the function of 
pulse. Here, the motion of the nature of the function of pulse may be a periodic motion in 
the general sense of the word. In this case, it is provided that the motion of the nature of 
the function of pulse produced in each cycle generated by the motion generating machine 
occurs within the period of such an extent as avoids exerting an influence on the 
measurement using the motion of the nature of the function of pulse produced in the next 
cycle. The aforementioned motion of the nature of the function of pulse may be followed 
by a motion for eliminating the displacement produced by the motion. 

The fourth aspect of the invention, besides the third aspect of the invention, is 
characterized by being provided with a first converting means for determining the Fourier 
component on the frequency axis of the motion of the nature of the function of pulse 
mentioned above and a second converting means for determining the Fourier component on 
the frequency axis of the output of an acceleration measuring device, an angular velocity 
measuring device, or an angular acceleration measuring device which is a device subject to 



calibration mentioned above and farther provided with a means to display or a means to 
transmit the frequency characteristic of the correction of an acceleration measuring device, 
an angular velocity measuring device, or an angular acceleration measuring device which is 
a device subject to calibration obtained from the outputs of the first and second converting 
means. 

The fifth aspect of the invention is characterized by the fact that the motion 
generating machine mentioned above is a motion generating machine which produces a 
random motion. The term "random motion" as used herein means that the motion can be 
handled as a white noise within the range of frequency band subject to calibration. 

The sixth aspect of the invention, besides the fourth aspect of the invention, is 
characterized by being provided with a first converting means for determining the Fourier 
component on the frequency axis of the random motion mentioned above and a second 
converting means for determining the Fourier component on the frequency axis of the 
output of an acceleration measuring device, an angular velocity measuring device, or an 
angular acceleration measuring device which is a device subject to calibration mentioned 
above and further provided with a means to display or a means to transmit the frequency 
characteristic of an acceleration measuring device, an angular velocity measuring device, or 
an angular acceleration measuring device which is a device subject to calibration obtained 
from the outputs of the first and second converting means. 

Before the measurement by the use of a laser interferometer is started, the control 
to the extent aimed at is possibly blocked as by the parasitic vibration of the motion 
generating machine. In this case, by feeding back the information showing the state of 
motion of the motion generating machine, it is made possible to attain the control as aimed 
at and cope with the trouble so as to repress the influence of the variation with time of the 
device. The seventh aspect of this invention, therefore, is characterized by the fact that 
the motion generating machine mentioned above is provided with an accelerometer for 
controlling a motion and a feedback circuit or a controlling device for controlling the drive 
device mentioned above so that the signal from the accelerometer may agree with a 
predetermined value. 

The eighth aspect of the invention is characterized by having as one object thereof 
the performance of the surveillance mentioned above with high precision and possessing a 
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structure for imparting dynamic matrix sensitivity to the aforementioned accelerometer for 
controlling a motion and, in the control of the feedback or in the control of a calculating 
machine, estimating the motion of the sample mounting table of the motion generating 
machine from the output vector of the accelerometer for controlling the motion by using 
the dynamic matrix sensitivity mentioned above, and controlling the motion. 

The ninth aspect of the invention is aimed at evaluating an error and is characterized 
by being provided with a means to display for determining an error from the output of the 
data processing device mentioned above and the output of an acceleration measuring 
device, an angular velocity measuring device, or an angular acceleration measuring device 
which is a device subject to calibration and displaying the value of this error or a means to 
convey for transmitting this value. 

Where the laser interferometer cannot be used easily, an alternate of high precision 
is used instead. This alternate at times abhors an appreciable degradation of measuring 
precision. The 10 th aspect of the invention is directed toward coping with this situation 
and is characterized by being provided with a motion generating machine capable of 
inducing a translational motion or a rotary motion, an acceleration measuring device, an 
angular velocity measuring device, or an angular acceleration measuring device which is a 
device subject to calibration and is at least temporarily fixed on the motion generating 
machine, an output means for taking out the output from the device subject to calibration, 
an inertia sensor capable of seizing a multidimensional motion calibrated in advance by 
using a dynamic matrix sensitivity measuring device of the inertia sensor recited in claim 1, 
a data processing device for processing the data showing the state of multidimensional 
motion obtained from the inertia sensor and converting a multidimensional translational 
motion or a multidimensional rotary motion to the value fit for the predetermined 
coordinate system, and a means to display or a means to convey the output of the data 
processing device and the output of the acceleration measuring device, the angular velocity 
measuring device, or the angular acceleration measuring device which is a device subject to 
calibration. 

When an acceleration measuring device, an angular velocity measuring device, or an 
angular acceleration measuring device is calibrated on the ground, since the calibration is 
carried out in the presence of gravitational acceleration, it is preferred to be capable of 
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excluding the influence thereof from the measured value. The 1 1 th aspect of the invention, 
therefore, is characterized by being provided with a means of adjusting the direction of 
mounting an acceleration measuring device, an angular velocity measuring device, or an 
angular acceleration measuring device which is a device subject to calibration so that the 
direction may be varied relative to the direction of gravity and a means for obtaining an 
output from the aforementioned device subject to calibration temporarily fixed as set in a 
plurality of directions relative to the direction of gravity, finding dynamic matrix sensitivity 
for each of the plurality of directions mentioned above, and estimating the dynamic matrix 
sensitivity liberated from the influence of gravity from the plurality of dynamic matrix 
sensitivities mentioned above. 

Particularly the aforementioned device subject to calibration which characterizes the 
12 th aspect of this invention is an acceleration measuring device, an angular velocity 
measuring device, or an angular acceleration measuring device which possesses the output 
thereof only in a specific direction. The motion generating machine which induces the 
aforementioned translational motion or rotary motion may be capable of producing a 
motion with a plurality of degrees of freedom. 

The 13 th aspect of the invention is characterized by the whole of a device for 
measuring the dynamic matrix sensitivity of an inertia sensor being particularly installed in a 
room capable of intercepting noise or vibration from the exterior, the device for measuring 
acceleration which is a device subject to calibration being a seismograph, and the motion 
generating machine being a motion generating machine which generates a motion of the 
nature of the function of pulse or the function of cycle and imparts small vibration in a 
vibration frequency band or in a seismic zone detectable by the seismograph (strong motion 
seismograph). The seismic wave is known in two kinds, a longitudinal wave and a 
transverse wave. Specifically, it is inherently a multidimensional motion which possibly 
occurs simultaneously in the vertical direction and the horizontal direction, which possibly 
occurs later in the horizontal direction, or which occurs in the form of a rotation of the 
ground and it is characterized by the fact that the direction of motion is an unknown factor. 
The motion mentioned above preferably conforms to this situation. 

The 14 th aspect of the invention is characterized by the aforementioned motion 
generating machine being particularly a motion generating machine which generates 
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random motion and imparts vibration in an vibration frequency band detectable by the 
inertia sensor used for constant detection of vibration of a vibro-isolating bed, and the 
motion generating machine or the vibro-isolating bed being provided with a cooling device 
or a thermostat. Here, the vibration of the floor surface for mounting the vibro-isolating 
bed comprises a horizontal motion and a vertical motion and is characterized by the fact 
that the direction of motion cannot be known in advance. Thus, the aforementioned 
motion-generating machine is preferred to be adapted to generate a motion conforming to 
this situation. 

The 1 5 th aspect of this invention is characterized by the aforementioned motion 
generating machine being particularly a motion generating machine which gives a vibration 
in a vibration frequency band detectable by an accelerometer used for controlling an 
automobile suspension or for controlling an automobile passenger protecting airbag, and 
the acceleration measuring device which is a device subject to calibration being an 
accelerometer used for controlling the automobile suspension or for controlling the 
automobile passenger protecting airbag and being provided with a temperature controlling 
device for controlling the temperature environment of the acceleration measuring device 
which is a device subject to calibration. Since the vibration of the automobile suspension 
is characterized by simultaneously producing a translational motion in the direction of a 
spring and a rotary motion around the rotational axis of the link of the suspension 
mechanism, the aforementioned motion generating machine is preferred to be enabled to 
generate a motion conforming to this situation. Further, the accelerometer to be used for 
controlling the automobile passenger protecting airbag has an important requirement of 
possessing the following characteristic features and the aforementioned motion generating 
machine is preferred to be capable of generating a motion conforming to this situation. 
First, the evaluation of safety laid out in the specification is stipulated to give results of 
calculation of the absolute values of acceleration in head-on collision and lateral collision 
below specified levels. The actual collision of cars does not always occur in the form of 
head-on collision. That is, the accelerometer which is used in the system for ensuring 
safety in collision must always detect correct acceleration without relying on the direction 
of acceleration The calibration, therefore, must be effected in terms of vector. 



The 16 th aspect of this invention is characterized by the aforementioned motion 
generating machine being particularly a motion generating machine which simultaneously 
generates a motion of the nature of the function of pulse or the nature of the function of 
cycle and comprising a translational motion and a rotary motion and imparts a motion in a 
vibration frequency band detectable by an inertia sensor used for controlling the motion of 
a robot and the aforementioned acceleration measuring device which is a device subject to 
calibration being an inertia sensor to be used for controlling the motion of the robot 
mentioned above. Since the motion of the robot is characterized by being capable of 
simultaneously generating a translational motion and a rotary motion with high accuracy, 
the motion mentioned above is preferred to conform to this situation. 

The 17 th aspect of this invention is characterized by the aforementioned motion 
generating machine being particularly a motion generating machine which gives a vibration 
in a vibration frequency band detectable by an inertia sensor to be used for measuring a 
human body motion, a vibration exerted on a human body, or an animal behavior monitor 
and the aforementioned acceleration measuring device which is a device subject to 
calibration being an inertia sensor, i.e. distributed accelerometers, to be used for measuring 
a human body motion, a vibration exerted on a human body, or an animal behavior monitor 
and being provided with a multichannel signal output terminal for emitting a signal for the 
distributed accelerometers. Since the vibration exerted on the inertia sensor to be used for 
measuring a human body motion, a vibration exerted on a human body, or an animal 
behavior monitor is characterized as indicated herein below, the aforementioned motion 
generating machine is preferred to be enabled to generate a motion conforming to this 
situation. 

1) The specification laid down by ISO regarding the measurement of vibration 
exerted on a human body (Non-Patent Document 6, Non-Patent Document 7 and 
Non-Patent Document 8) defines use of a three-axis accelerometer for vibrations in hands 
and arms and use of a six-axis accelerometer for whole body vibrations. To be specific, 
since the process of judgment embraces an operation of exerting a weight on the X 
component, Y component, and Z component of given acceleration, squaring the respective 
responses, totaling the produced squares, and reducing the total to a square root (raising to 
the l/2nd power), the specification evidently expects the measurement of the acceleration 
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in terms of vector. The value of weight differs between whole body vibration and 
vibration of hands and arms. Meanwhile, the method for calibrating an accelerometer 
concerns the calibration of the amplitude of a one-axis vibration. Evidently, the human 
body possesses an ability to discern a direction of vibration. The fact that the 
accelerometer which is used for investigating the influence of the vibration on a human 
body is calibrated exclusively with the amplitude of acceleration is irrational. 

2) The measurement of the motion of a human body yields a valuable basic data for 
the sports engineering and for the control of a humanoid robot. It is evidently irrational, 
however, to assume that the direction of a motion of a human body is constantly fixed. 
The fact that the accelerometer which is used for investigating the motion of a human body 
is calibrated exclusively with the amplitude of acceleration is irrational. 

3) The animal inhabits a real space and the translational motion occurs 
three-dimensionally. When the rotational motion is also taken into account, it is logical to 
conclude that the animal inhabits a six-dimensional space. The fact that the inertia sensor 
for monitoring a behavior is not calibrated in six-dimensional space may well be judged 
senseless. 

4) In the measurement of the motion of a human body and the analysis of a human 
body motion in the sport engineering, while inertia sensors are attached to joints, arms, and 
feet, the calibration must be made based on the motion of six degrees of freedom because 
the motion simultaneously induces a vertical and a horizontal motion and entails a rotary 
motion produced in a joint. 

The 18 th aspect of the invention, besides the 10 th aspect of the invention, is 
characterized by particularly the motion generating machine which induces the 
aforementioned translational motion or rotary motion generating a motion along one axis, 
the acceleration measuring device, the angular velocity measuring device, or the angular 
acceleration measuring device which is a temporarily fixed device subject to calibration 
being furnished with a one-axis output, and the output from the device subject to 
calibration being furnished with a multi-axis output. 

With the exception of the case in which the aforementioned device subject to 
calibration is one-axis, the accelerometer is expected to be calibrated with a greater degree 
of freedom than the degree of freedom possessed by the motion generating machine. The 
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number of axes which can be handled is possibly smaller than the number of output axes of 
the accelerometer expected to be used for calibration. The 19 th aspect of the invention, 
which concerns a method to be applied in such a case to the device for measuring the 
dynamic matrix sensitivity of the inertia sensor possessing the aforementioned 
characteristics, is characterized, on the assumption that the aforementioned displacement 
measuring means produces outputs of different N axes, the degree of freedom of the 
motion of the aforementioned motion generating machine is a natural number M of not 
more than 6, and the relation M x N = 1 is not satisfied, by dividing the M degrees of 
freedom into a plurality of groups allowed duplication, 

1) obtaining input vectors and corresponding output vectors for each of the groups 
mentioned above by measurement, 

2) subjecting the input vectors and the output vectors to conversion and integration 
into input vectors possessing M degrees of freedom and corresponding output vectors, and 
subsequently 

3) deriving a dynamic matrix sensitivity from the correspondence of the integrated 
output vectors to the integrated input vectors. 

As other characteristics of the invention, particularly the aforementioned motion 
generating machine is a motion generating machine which gives vibration in a vibration 
frequency band detectable by an inertia sensor to be used in an inertial navigation system 
mounted on an automobile, a submarine, a torpedo, a missile, an airplane, or a guiding 
device therefor and the aforementioned acceleration measuring device which is a device 
subject to calibration is an inertia sensor to be used in an inertial navigation system 
mounted on an automobile, a submarine, a torpedo, a missile, an airplane, or a guiding 
device therefor. Since the vibration which is exerted on the inertia sensor to be used in an 
inertial navigation system mounted on an automobile, a submarine, a torpedo, a missile, an 
airplane, or a guiding device therefor is characterized by simultaneously giving rise to a 
horizontal motion, a vertical motion, and a rotary motion, the aforementioned motion 
generating machine is preferred to be enabled to generate a motion conforming to these 
characteristics. 

As other characteristics of the invention, the aforementioned motion generating 
machine is particularly a motion generating machine which gives vibration in a vibration 
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frequency band detectable by an inertia sensor to be used for preventing the image of an 
image pickup device or a motion picture device from blurring and the aforementioned 
acceleration measuring device which is a device subject to calibration is an inertia sensor to 
be used for preventing the image of an image pickup device or a motion picture device 
from blurring. Here, since the prevention of the image of an image pickup device or a 
motion picture device from blurring has an important requirement of fulfilling simultaneous 
detection of a horizontal motion, a vertical motion, and a rotary motion, the 
aforementioned motion generating machine is preferred to generate a motion conforming to 
this situation. 

As other characteristics of the invention, the aforementioned motion generating 
machine is particularly a motion generating machine which gives vibration in a vibration 
frequency band detectable by an accelerometer to be used for controlling an automobile 
passenger protecting airbag and the aforementioned acceleration measuring device which is 
a device subject to calibration is an accelerometer to be used for controlling the automobile 
passenger protecting airbag. Here, since the accelerometer to be used for controlling the 
automobile passenger protecting airbag has important characteristic features of fulfilling the 
following requirement, the aforementioned motion generating device is preferred to be 
enabled to generate a motion conforming to this situation. 

First, though the results of calculation of the absolute values of acceleration in the 
head-on collision and the lateral collision obtained in the evaluation of safety stipulated by 
the specification are expected to be below a prescribed level, the actual collisions of cars do 
not always occur in the form of head-on collision or lateral collision. The accelerometer 
to be used in the collision safety system, therefore, must always detect the acceleration 
correctly without relying on the direction of exertion of the acceleration. That is, the 
calibration must be made in terms of vector. 

On the assumption that the inertia sensor possesses K axes, the output vector 
possibly form a K-dimensional space, and the motion generating machine possesses M 
degrees of freedom, the K x M-dimensional matrix can be derived for the inertia sensor as a 
sensitivity matrix. The act of solving all the components of the sensitivity matrix as the 
function of frequency or as the function of the frequency and the environmental parameters 
such as the temperature of the inertia sensor and the posture relative to the gravitational 
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acceleration shall be regarded as the calibration. 

Since all the components of this matrix are unknown and since the motion 
generating machine having the inertia sensor mounted thereon is capable of generating a 
motion of M degrees of freedom, by inducing generation of M independent vector motions 
and measuring the produced output signals from the K axes, it is made possible to 
formulate simultaneous linear equations concerning all the components of the K x 
M-dimensional sensitivity matrix. The solution of the simultaneous linear equations 
results in deriving a sensitivity matrix composed of the normal sensitivities and the cross 
sensitivities in a given angular frequency co. 

Even by using a one-axis motion generating machine (actuator), it is theoretically 
possible to form linearly independent excitation vectors as viewed from the inertia sensor 
which is subject to calibration. This formation, however, necessitates a variation in the 
attachment of the multi-axis inertia sensor. The care should be taken to the fact that the 
detecting mechanisms involved in a plurality of detecting axes do not necessarily function 
simultaneously while they constantly maintain the influence of gravitational acceleration 
uniformly without being affected by the attachment with jigs. 

Brief Description of the Drawings: 

Fig. 1 is a type section showing the configuration for embodying this invention. 

Fig. 2 is a diagram showing the acceleration exerted on the conventional 
accelerometer and the resolution of the vector thereof 

Fig. 3 is a side view illustrating an accelerometer set on the conventional one-axis 
motion generating machine. 

Fig. 4 is a side view illustrating a conventional reference accelerometer and an 
accelerometer subject to calibration, which are joined in series connection with the 
sensitivity axis aligned with the direction of motion generated by a motion generating 
machine. 

Fig. 5 is a diagram showing an assembly of one-axis inertia sensors and the relation 
between the matrix sensitivity of the assembled inertia sensors and the matrix sensitivity of 
each of the inertia sensors. 
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Fig. 6 is a perspective view showing the coordinate system fixed on an 
accelerometer subject to the calibration contemplated by this invention. 

Fig. 7 is a perspective view illustrating the coincidence between the coordinate axis 
ZX of a motion generating machine and the ZaXa axis of an accelerometer according to 
this invention. 

Fig. 8 is a perspective view illustrating the coincidence between the coordinate axis 
ZX of a motion generating machine and the ZaYa axis of an accelerometer according to 
this invention. 

Fig. 9 is a type section illustrating the configuration for embodying this invention. 
The reference numerals in the diagrams denote the following components. 

1 Accelerometer subject to calibration 

2 Controlling or measuring real time laser interferometer (laser) 

3 Controlling or measuring real time laser interferometer (corner cube) 

4 Controlling inertia sensor having the matrix sensitivity defined 

5 Horizontal translational motion actuator and rotary motion generating 
actuator 

6 Vertical translational motion actuator and rotary motion generating actuator 

7 Horizontal translational motion actuator and rotary motion generating 
actuator 

8 Power amplifier 

9 Processor (CPU) 

10 Cooling water piping system 

11 Table 

Best Mode for carrying out the Invention: 

This invention will be described more specifically below based on working examples 
thereof with reference to the accompanying drawings. 

Generally, on the assumption that an inertia sensor possesses K sensitivity axes and 
an output vector space possibly possesses K dimensions and that a motion generating 
machine as a source of vibration possesses M degrees of freedom, a K x M-dimensional 
matrix can be derived as the sensitivity matrix of the inertia sensor. The act of solving all 
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the components of the sensitivity matrix as the function of frequency or as the function of 
the frequency and the environmental parameters such as the posture and the temperature of 
the inertia sensor relative to the gravitational acceleration is regarded as calibration. 

Since all the components of this matrix are unknown and since the motion 
generating machine having the inertia sensor mounted thereon is capable of generating a 
motion of M degrees of freedom, by inducing generation of M independent vector motions 
and measuring the produced output signals from the K axes, it is made possible to 
formulate simultaneous linear equations concerning all the components of the K x 
M-dimensional matrix. The solution of the simultaneous linear equations results in 
deriving a sensitivity matrix composed of the normal sensitivities and the cross sensitivities 
in a given angular frequency a>. 

Now, the components of one example of the configuration for embodying this 
invention will be explained below with reference to Fig. 1 . 

Fig. 1 illustrates an accelerometer 1 as a device subject to calibration, a laser 2 for a 
controlling or measuring real time laser interferometer, a corner cube 3 for the controlling 
or measuring real time laser interferometer, an inertia sensor 4 for controlling a motion 
generated by a motion generating machine, a horizontal translational motion actuator and a 
rotational motion generating actuator 5 and 7 as sources of vibration, a vertical 
translational motion actuator and rotary motion generating actuator 6 as a source of 
vibration, a power amplifier 8, a processor 9 concurrently serving as a data processing 
device, a cooling water piping system 10, and a table 11 as the base for attachment of an 
inertia sensor subject to calibration. 

Here, the laser 2 of the controlling or measuring real time laser interferometer and 
the corner cube 3 which constantly returns the laser beam from the laser 2 in parallel with 
the incident light integrally form a controlling or measuring real time laser interferometer 
and the inertia sensor 4 for controlling a motion generated by a motion generating machine 
has a matrix sensitivity defined therein. As the device subject to calibration, an angular 
velocity measuring device or an angular acceleration measuring device may be selected 
besides the aforementioned accelerometer as an acceleration measuring device. 

Further, this device subject to calibration does not need to be a simple acceleration 
measuring device, angular velocity measuring device, or angular acceleration measuring 
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device but may be a composite body formed by combining such devices. For example, the 
composite body may be formed by having one-axis accelerometers disposed severally along 
the X, Y and Z axes of the cubic coordinates or the R, 0 and q> axes of the spherical 
coordinates. Fig. 5 illustrates an example of the accelerometer having one-axis 
accelerometers disposed severally along the X, Y and Z axes. 

The configuration of Fig. 1 uses a laser interferometer. The fact that a reference 
accelerometer calibrated with a laser interferometer can be used to calibrate other 
accelerometers in place of the laser interferometer is evident from the layout and 
configuration of Fig. 1. 

Now, the best mode of embodiment of this invention will be explained below based 
on working examples. 
(Example 1) 

First, as Example 1, the case of adopting a one-axis accelerometer as the 
accelerometer 1 subject to calibration and using a three-dimensional motion generating 
machine will be explained. 

The operating part of the processor 9 is manipulated to transmit control signals to 
the individual actuators 5, 6 and 7 which are sources of vibration so that the tables of these 
motion generating machines may generate respectively prescribed motions. At this time, 
the inertia sensor 4 for controlling the motion generating machine makes measurement 
through the table 1 1 to judge whether or not the individual actuators 5, 6 and 7 are 
operating with prescribed vibrations and the processor 9 generates a control signal based 
on the difference between the waveform of the target motion of the table and the actually 
measured waveform and transmits control signals to the individual actuators 5, 6 and 7 so 
that the actuators 5, 6 and 7 may produce prescribed vibrations. 

Thus, the prescribed vibrations are imparted to the actuators 5, 6 and 7 and the 
output signal from the accelerometer 1 as an inertia sensor is input to the processor 9 by 
using one of the sensitivity axes of the accelerometer 1 as the X axis. Here, for the sake 
of brevity, the signals produced by carrying out Laplace transform to the output signals will 
be used in the following explanation. The output signals which have undergone the 
Laplace transform will be denoted by (a^co), 0, 0). The letter co denotes the angular 
frequency. The Y axis component and the Z axis component are both zero as a matter of 
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course. 

The controlling or measuring real time laser interferometers 2 and 3, which are 
means for measuring the input to the accelerometer 1, measure the input acceleration to the 
accelerometer 1 through the table 11. The measured input acceleration which has 
undergone the Laplace transform are denoted by (a^QcD), &i y Q(o\ a^Qo))). It is not 
stipulated that the input acceleration falls on the sensitivity axis of the accelerometer 1. 
When this stipulation is accepted, it does not result in supporting a conclusion that the 
input acceleration is regarded as a vector. 

Incidentally, the Laplace transform of the input acceleration is used here for the 
purpose of specifically indicating the response to the input in accordance with the definition 
of transfer function. 

Thus, the actuators 5, 6 and 7 are exposed to the respective prescribed vibrations 
and the output signals from the accelerometer 1 are input to the processor 9 by using one 
of the sensitivity axes of the accelerometer 1 as an inertia sensor as the X axis. Here, the 
output signals which have undergone the Laplace transform according to the definition of 
transfer function will be used in the explanation. The output signals which have 
undergone the Laplace transform will be denoted by (aoxQ&X 0> 0). The letter co denotes 
the angular frequencies. The Y axis component and the Z axis component are both zero 
as a matter of course. 

The controlling or measuring real time laser interferometers 2 and 3 that are means 
for measuring the input to the accelerometer 1, measure the input acceleration to the 
accelerometer 1 through the table 11. The measured input acceleration which has 
undergone the Laplace transform are denoted by (aixGa>)> a^Qo)), Mjg))). It is not 
stipulated that the input acceleration falls on the sensitivity axis of the accelerometer 1. 
When this stipulation is accepted, it does not result in supporting a conclusion that the 
input acceleration is regarded as a vector. 

Incidentally, the Laplace transform of the input acceleration is used here for the 
purpose of specifically indicating the response to the input in accordance with the definition 
of transfer function. The method for numerically processing the acquired data is known in 
various kinds and the software for this processing is easily procurable. Thus, the 
numerical data aimed at may be determined based on the condition of the data and the 
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target precision of processing. 

In this case, the matrix sensitivity is expressed by the matrix of 1 x 3 (Sx^, S^y, Sxz). 
Sx^c, Sx, y and Sxz denote respectively the normal sensitivity expressing the relation of the X 
axis output signal to the X axis input component of the accelerometer 1, the cross 
sensitivity expressing the relation of the X axis output signal to the Y axis input component 
of the acceleration, and the cross sensitivity expressing the relation of the X axis output 
signal to the Z axis input component of the acceleration. At this time, the relation 
between the output signal and the input signal is expressed by the aforementioned 
Mathematical 1. 

The three-dimensional motion generating machine is capable of forming three 
independent vector motion accelerations in a three-dimensional space. These vector 
motion accelerations are denoted respectively by (a^i, ai y i, aizi)exp(jcot), (aix2, ai y2 , 
ai Z 2)exp(joot) and (aix3, ai y3j a^expG'cDt). By substituting these accelerations in the 
Mathematical 1, the following relation is established. The terms aoxi, aoyi, and aozi denote 
complex numbers as the functions of co taking into account the initial phase relative to 
exp(joot) and the gain. 

(Mathematical 2) 

«atl(^)=5 F x.x(^)«ft1 +Sx,y(*>y*iyi 

a cot2( a} ) = S x f x( a) ) a ix2 + S x,y( <D ) a iy2 +S x,z( m ) a ti2 

a cx3( a) ) = SxjA* )*r*3 + S x % y{ a )°y3 +S x,z( a, ) a ii3 

Mathematical 2 represents simultaneous linear equations which concern (S^, Sx,, y , 
Sx^) in each of the numbers of angular frequency co and the coefficient matrix has a solution 
with respect to (S^ Sx,, y , Sx^) because (aixk, ai yk , a^) (k = 1, 2, 3) are independent and, 
consequently, the determinant is not zero. 

The consequence is that the normal sensitivity Sx^, and cross sensitivities Sx, y and 
Sx^ are determinable. By acquiring an experimental data with various numbers of angular 
frequencies and solving the simultaneous linear equations of Mathematical 2 obtained with 
the various angular frequencies, it is made possible to measure the normal sensitivity and 
the lateral sensitivity as the functions of the angular frequencies. 
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(Example 2) 

Example 2 represents the case of adopting a two-axis accelerometer as the 
accelerometer 1 subject to calibration and using a three-dimensional motion generating 
machine. The output signals of the accelerometer 1 are input to the processor 9 by using 
the two sensitivity axes of the accelerometer 1 as X axis and Y axis. The processor 9 
subjects the output signals to Laplace transform and consequently obtains (ao X (ja>), e^yQco), 
0). The letter go denotes the angular frequencies. The Z axis component is zero of 
course. Similarly, the laser interferometer 2 measures the input acceleration of the 
accelerometer 1 and the processor 9 subjects the input acceleration to Laplace transform 
and consequently obtains (a^Oco), ai y (jco), aiz(jco)). It is not stipulated that the input 
acceleration falls on a plane (sensitivity plane) determined by the two sensitivity axes of the 
accelerometer 1. When this stipulation is accepted, it does not result in supporting a 
conclusion that the input acceleration is regarded as a vector. At this time, the matrix 
sensitivity is expressed by the matrix of 2 x 3 shown by the following Mathematical 3. 

(Mathematical 3) 



S,^, Sx, y and Sx, 2 denote respectively the normal sensitivity expressing the relation of 
the X axis output signal to the X axis input component of the accelerometer 1, the cross 
sensitivity expressing the relation of the X axis output signal to the Y axis input component 
of the acceleration, and the cross sensitivity expressing the relation of the X axis output 
signal to the Z axis input component of the acceleration. Each of S y>x , S y , y and S y ^ 
expresses the relation of the input signal to the Y axis output of the accelerometer 1. S y>x 
represents the cross sensitivity expressing the relation of the Y axis output signal to the X 
axis input component of the acceleration, S y , y represents the normal sensitivity expressing 
the relation of the Y axis output signal to the Y axis input component of the acceleration, 
and S y , z represents the cross sensitivity expressing the relation of the Y axis output signal to 
the Z axis input component of the acceleration. In this case, the relation between the 
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input vector and the output vector is expressed by Mathematical 4. 
(Mathematical 4) 



f a m {a>) exp(jat) > _ ( S xx S^ y (a>) , z (co) 
a<p (co) exp(Jo t))[s ytX (o\ (co), S^ 2 (co) 



OfyByqpijcDt) 
X^&xpijcot) j 



The three-dimensional motion generating machine is capable of forming three 
independent vector vibration accelerations in a three-dimensional space. These 
accelerations will be expressed as follows: 

(aixi, ai y i, aizi)exp(j(Dt), 

(aix2, ai y2 , aiz2)exp(jo)t) and 

(aixs, ai y3 , a^expGcat). 

By substituting these accelerations in Mathematical 1, the following relation is 
established. The terms ao X i and ao y i (i = 1, 2, 3) are complex numbers as the functions of co 
which take into account the initial phase relative to exp(jcot) and the gain. When these 
complex numbers are substituted in Mathematical 4, the following three formulas, i.e. 
Mathematical 5, Mathematical 6 and Mathematical 7, are established. When these 
formulas are unified as simultaneous linear equations concerning all the components of the 
sensitivity matrix, Mathematical 8 is derived. 

(Mathematical 5) 
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(Mathematical 6) 



' *a*2(*0 expO© f) W S^a), S Xt y(a>),S XiZ (a>) ^ 



(Mathematical 7) 



^expOcy*)' 



(Mathematical 8) 
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Since the three vectors (a^, a iyk , a^) (k = 1, 2, 3) in Mathematical 8 are linearly 
independent, it is only plain that the six vectors shown below: 

(aixk, a iyk , ai*, 0, 0, 0) (k = 1, 2, 3) and 

(0, 0, 0, aixk, ^ a^) (k = 1 ? 2, 3) 
are linearly independent. Since the determinant of the simultaneous linear equations is not 
zero, it is only logical to conclude that they are capable of being solved without fail. The 
sensitivity Sij (i = x, y, j = x, y, z) is a normal sensitivity when the suffixes are equal or a 
cross sensitivity when the suffices are not equal. The processor 9 is enabled, by acquiring 
an experimental data with a varying angular frequencies and solving eight sets of 
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simultaneous linear equations obtained individually with the varying numbers of angular 
frequencies, to measure the normal sensitivities and the cross sensitivities as functions of 
the angular frequencies. 
(Example 3) 

Example 3 represents the case of adopting a three-axis accelerometer as the 
accelerometer 1 (inertia sensor) subject to calibration and using a three-dimensional motion 
generating machine. The output signals of the accelerometer 1 are input to the processor 
9 by using the three sensitivity axes of the accelerometer 1 as the inertia sensor, i.e. X axis, 
Y axis and Z axis. The processor 9 subjects the output signals to Laplace transform and 
consequently obtains (&o K (}(o\ aoy(ja>), aozG* 0 ))- The letter go denotes the angular 
frequencies. 

Similarly, the laser interferometer 2 measures the input acceleration of the 
accelerometer 1 and the processor 9 subjects the input acceleration to Laplace transform 
and consequently obtains (aix(ja>), a^co), a^co)). It is not stipulated generally that the 
input acceleration falls on a vector space (sensitivity space) determined by the three 
sensitivity axes of the accelerometer. When this stipulation is accepted, it does not result 
in supporting a conclusion that the input acceleration is regarded as a vector. The 
accelerometer 1 projects the input acceleration space on a sensitivity vector space. At this 
time, the matrix sensitivity is expressed by the matrix of 3 x 3 shown below. 

(Mathematical 9) 



^ e o c ^ 
^x.x^x.y^x.z 

e <? 



J 

Each of Sx^c, Sx, y and Sx, z expresses the relation of the input signal to the 
accelerometer 1 to the X axis output of the accelerometer 1 . Sx^ denotes the normal 
sensitivity expressing the relation of the X axis output to the X axis input component, Sx, y 
denotes the cross sensitivity expressing the relation of the X axis output signal to the Y axis 
input component of the acceleration, and Sx, z denotes the cross sensitivity expressing the 
relation of the X axis output signal to the Z axis input component of the acceleration. 
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Each of S ytX , Sy, y and S y ^ expresses the relation of the input signal to the Y axis output of 
the accelerometer 1. S y ^ denotes the cross sensitivity expressing the relation of the Y axis 
output signal to the X axis input component of the acceleration, S y , y denotes the normal 
sensitivity expressing the relation of the Y axis output signal to the Y axis input component 
of the acceleration, and S y , z denotes the cross sensitivity expressing the relation of the Y 
axis output signal to the Z axis input component of the acceleration. Each of S^, Sz,, y and 
Sz^ expresses the relation of the input-output signal to the Z axis output of the 
accelerometer 1. 

Sz^ denotes the cross sensitivity expressing the relation of the Z axis output signal 
to the X axis input component of the accelerometer 1, Sz, ?y denotes the cross sensitivity 
expressing the relation of the Z axis output signal to the Y axis input component of the 
accelerometer 1, and Sz^ denotes the normal sensitivity expressing the relation of the Z axis 
output signal to the Z axis input component of the accelerometer 1. At this time, the 
relation between the input vector and the output vector is expressed by the following 
Mathematical 10. 

(Mathematical 10) 



a oy (a>)exp(jat) 
a m Uo)&yp{ja>t)j 



{S zx {co),S zy {a>),S zz {a>) j 



V 



The three-dimensional motion generating machine is capable of forming three 
independent vector vibration accelerations in a three-dimensional space. These 
accelerations will be expressed as (aixi, a iy i, aizOexpQcot), (a^, ai y2 , a^expQcot), and (a^, 
ai y3 , a i23 )exp(j(Dt). 

By substituting these accelerations in Mathematical 10, the following three 
equations, Mathematical 11, Mathematical 12 and Mathematical 13 are established. The 
terms a^acy* and a^i = 1, 2, 3) are complex numbers as the functions of a> which take 
into account the initial phase relative to expQcot) and the gain. When these formulas are 
unified as simultaneous linear equations concerning all the components of the sensitivity 
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matrix, Mathematical 14 is derived. 
(Mathematical 11) 



a m l(a>)expOatf) 



(Mathematical 12) 



'a m2 {a>)«up{jat) y 



(Mathematical 13) 



S x ,x(<l>hSx,yW.Sx.z (*>) 
[S z x (c0),S zy (a>),S Z 2 (a>) } 



ct iy ^exp(ja>t) 



' S XtX {a>),S xy {&),S xa {(t>) ^ 
K S ZiX {a>),S Zty (co),S Z 2 (a>) j 



^2 exp(ja)t)j 



(Mathematical 14) 



[S zx (a>),S zy (a>),S zz (a>) J 
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Since the three vectors (a^, ai yk , a^) (k = 1, 2, 3) in Mathematical 14 are linearly 
independent, it is only plain that the 9 vectors shown in the coefficient determinant 
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(aixk, a iyk , aizk , 0, 0, 0, 0, 0, 0) (k = 1, 2, 3) 

(0, 0, 0, a^, a iyk , 0, 0, 0) (k = 1, 2, 3) and 

(0 5 0, 0, 0, 0, 0, aixk, ai^, %*) (k = 1, 2, 3) 

are linearly independent. Since the coefficient determinant of Mathematical 14 is not zero, 
it is only logical to conclude that they are capable of being solved without fail. The 
sensitivity Sg (i = x, y, z, j = x, y, z) is a normal sensitivity when the suffixes are equal or a 
cross sensitivity when the suffices are not equal. The processor 9 is enabled, by acquiring 
an experimental data with a varying angular frequencies and solving the simultaneous linear 
equations of Mathematical 14 obtained individually with the varying angular frequencies, to 
measure the normal sensitivities and the cross sensitivities as functions of the angular 
frequencies. 
(Example 4) 

Example 4 represents the case of adopting a three-axis accelerometer + an angular 
accelerometer of one degree of freedom as an accelerometer 1 (inertia sensor) subject to 
calibration and using a four-dimensional motion generating machine (translational motion 
of 3 degrees of freedom and angular acceleration of one degree of freedom). The 
sensitivity axes of the three translational motions of the accelerometer 1 as an inertia sensor 
are denoted by X axis, Y axis and Z axis and the input axis and the output axis of the 
angular acceleration around the X axis are denoted by a axis. 

The output signals of the accelerometer 1 are input to the processor 9 and the 
processor 9 subjects the output signals to Laplace transform and consequently obtains 
(a<*G<D), aoy(ja>), ao Z (jco) and aoaQ®))- The letter a> denotes the angular frequencies. 
Similarly, the laser interferometer 2 measures the input acceleration of the accelerometer 1 
and the processor 9 subjects the input acceleration to Laplace transform and consequently 
obtains (a^Qcc), ai y (ja>), a^co), ajaQo)). 

It is not stipulated generally that the input acceleration falls on a vector space 
(sensitivity space) determined by the four sensitivity axes of the inertia sensor. When this 
stipulation is accepted, it does not result in supporting a conclusion that the input 
acceleration is regarded as a vector. The accelerometer 1 is a sensor which fulfills the 
function of projecting the input acceleration space on the output acceleration space. At 
this time, the matrix sensitivity is expressed by the matrix of 4 x 4 as shown below. 



-28- 



(Mathematical 15) 

* S y*y 5lS >.* ? ^>.« 



Each of Sx, x , Sx, y , Sx, z and S^ a expresses the relation of the input signal to the X 
axis output of the accelerometer 1 . Sx, x denotes the normal sensitivity expressing the 
relation of the X axis output to the X axis input component, Sx, y denotes the cross 
sensitivity expressing the relation of the X axis output signal to the Y axis input component 
of the acceleration, Sx, z denotes the cross sensitivity expressing the relation of the X axis 
output signal to the Z axis input component of the acceleration, and Sx, a denotes the cross 
sensitivity expressing the relation of the X axis output signal to the a axis input component 
of the acceleration. 

Each of S y ,x, Sy, y , S y ^ and S y , a expresses the relation of the input signal to the Y 
axis output of the accelerometer 1. S y>x denotes the cross sensitivity expressing the 
relation of the Y axis output signal to the X axis input component of the acceleration, S y , y 
denotes the normal sensitivity expressing the relation of the Y axis output signal to the Y 
axis input component of the acceleration, S y ^ denotes the cross sensitivity expressing the 
relation of the Y axis output signal to the Z axis input component of the acceleration, and 
Sy, tt denotes the cross sensitivity expressing the relation of the Y axis output signal to the a 
axis input component of the acceleration. 

Each of Sz^c, Sz,y 5 Sz^ and Sz, a expresses the relation of the input signal to the Z axis 
output of the accelerometer 1 . Sz, x denotes the cross sensitivity expressing the relation of 
the Z axis output signal to the X axis input component of the acceleration, S^y denotes the 
cross sensitivity expressing the relation of the Z axis output signal to the Y axis input 
component of the acceleration, Sz, z denotes the normal sensitivity expressing the relation of 
the Z axis output signal to the Z axis input component of the acceleration, and Sz, a denotes 
the cross sensitivity expressing the relation of the Z axis output signal to the a axis input 
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component of the acceleration. 

Each of Scvc, So, y , So^ and So,a expresses the relation of the input signal to the a axis 
output of the accelerometer 1 . So, x denotes the cross sensitivity expressing the relation of 
the a axis output signal to the X axis input component of the acceleration, So, y denotes the 
cross sensitivity expressing the relation of the a axis output signal to the Y axis input 
component of the acceleration, Sa^ denotes the cross sensitivity expressing the relation of 
the a axis output signal to the Z axis input component of the acceleration, and So, a denotes 
the normal sensitivity expressing the relation of the a axis output signal to the a axis input 
component of the acceleration. At this time, the relation between the input vector and the 
output vector is expressed by the following Mathematical 16. 

(Mathematical 16) 

The four-dimensional motion generating machine is capable of forming four 
independent vector vibration accelerations in a four-dimensional space. These 
accelerations will be expressed as follows: 

(aixi, %i, ajzi, aiai)exp(jcot), 

(aM, ai y2 , a^, aia2)exp(jcot), 

(aixs, %3, a^, a ia 3)exp0cot) and 

(aix4 5 ai y4 , a^, aia4)exp(ja>t). 

By substituting these accelerations in Mathematical 16, the following four equations, 
Mathematical 17, Mathematical 18, Mathematical 19 and Mathematical 20 are established. 
The terms a 0X i,a 0 yj, aozi and a<>ai(i = 1, 2, 3, 4) are complex numbers as the functions of co 
which take into account the initial phase relative to exp(ja)t) and the gain. When the four 
equations, i.e. Mathematical 17, Mathematical 18, Mathematical 19 and Mathematical 20, 
are unified as simultaneous linear equations concerning all the components of the sensitivity 
matrix, Mathematical 21 is derived. 
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(Mathematical 17) 

^(©JexpO©* ) ^ ;r (a)),5 v> ,(®),5' > , r (a)),5 > ^( < 2>) ^expCfarf) 

,*«*l<<»)expO<»* )J [^^(^^^(^.^(^.^(^jJUfaiex^Oj 

(Mathematical 18) 

'a m2 (<D)exp<ija>0^ (S„(a>) S xy (,co) S XtI (a>) S,, a (tf>)Y a U2 exp<jc»0^ 

a^ 2 (o)expO^O = ^(o) S,,:(°>) S ya (eo) o^expO©/) 

a (a2 (G))expO < »0 5^(a>) S,j(a>) S ZrZ (co) S za (a>) a^expUat) 

< « < , a2 (^)expOo/)J [SaA^ S ay (a>) S az (co) 5 aia (fi?)J^a iff2 expC/a>/), 

(Mathematical 19) 

' a mZ {e>)BXv(ja>t)\ f^^X^^).^^)^^®) Va^expOof)^ 

a «3( a, > ex P0 fl, 0 5 *^>'^>»>. 5 * > x< fl> >.'W fi, > a^Wtpijot) 



(Mathematical 20) 

v {a>)exp(>0 ^(o), ^^(fi)),^^ z (fi»),5 >ra (fi>) a ^expO»/) 
^4(»)«p0»0 ~ 5 ><ir ((»> l 5 aof (<»),^ v (4»> l 5 r<a (d») < 7jfe4 expO»0 

^«*4(<»>«CPO<»OJ [ff«*(»>>*kr<">t*<U (®). W°>> JUto4 «PlA»0, 



(Mathematical 21) 
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' S, Si s»i 0 0 0 0 0 0 0 0 0 0 0 0 

0 0 0 0 Si Si s«i 0 0 0 0 0 0 0 0 

0 0 0 0 0 0 0 D ^ ^ ^ ^ 0 o 0 0 

0 0 0 0 0 0 0 0 0 0 0 0 Si Si a *i ^roJ 

^2 Ss ^ 3 a ia2 oooooooooooo 

0 0 0 0 s* S2 «2f«2 ^oa 0 0 0 0 0 0 0 0 

0 0 0 0 0 0 0 0 (3^2 Sa <2«2 S*2 0 0 0 0 

0 0 0 0 0 0 0 0 0 0 0 0 <2» 2 S2 «*2 ^2 

^3 S3 ^3 <^e3 0 oo 0 0 0 0 0 0 0 0 0 

0 0 0 0 <z iT3 S3 «2«3 0 0 0 0 0 0 0 0 

0 0 0 0 0 0 0 0 S3 S3 S 3 S*3 o 0 0 0 

0 0 0 0 0 0 0 0 0 0 0 0 S3 S3 ^3 S*3 

S4 s* a f*4 oooooooooooo 

o o o o sa s* ^oa oooooooo 

0 0 0 0 0 0 0 0 s 4 S4 *m <2,«4 o 0 0 0 

s 0 0 0 0 0 0 0 0 0 0 0 0 s< <h,A S4 Sa4. 




Since the four vectors (a^, a iy k 5 a^*, a ia k) (k = 1, 2, 3, 4) in Mathematical 21 are 
linearly independent, it is only plain that the 16 vectors shown in the coefficient determinant 
(aixk, a^, a^, a*,*, 0, 0, 0, 0, 0 5 0, 0, 0, 0, 0, 0, 0) (k = 1, 2, 3, 4), 
(0, 0, 0, 0, a^, aiyk, ai2k , a iak , 0, 0, 0, 0, 0, 0, 0, 0) (k = 1, 2, 3, 4), 
(0, 0, 0, 0, 0, 0, 0, 0, a^, a iyk , a^, a^, 0, 0, 0, 0) (k = 1, 2, 3, 4) and 
(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, aixk, a^, a^, a^) (k = 1, 2, 3, 4) 
are linearly independent. Since the coefficient determinant of Mathematical 21 is not zero, 
it is only logical to conclude that they are capable of being solved without fail. 

The sensitivity Sij (i = x, y, z, a, j = x, y, z, a) is a normal sensitivity when the 
suffixes are equal or a cross sensitivity when the suffices are not equal. This means that 
the normal sensitivity Sx^ and the cross sensitivities Sx, y , S*^ and S x « are obtained. The 
processor 9 is enabled, by acquiring an experimental data with a varying angular 
frequencies and solving the simultaneous linear equations of Mathematical 21 obtained 
individually with the varying numbers of angular frequencies, to measure the normal 
sensitivities and the cross sensitivities as functions of the angular frequencies. 
(Example 5) 

Finally, Example 5 adopts a three-axis accelerometer as an accelerometer 1 (inertia 
sensor) subject to calibration and uses a two-dimensional motion generating machine 
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(translation^ motion of two degrees of freedom) as shown in examples of the layout shown 
in Fig. 6 - Fig. 8. These diagrams illustrate the setups for embodying a method of 
calibration for the sensitivity of 3 x 3 matrix sought in the same manner as when an inertia 
sensor has three axes, and the motion vector space by the motion generating machine has 
three dimensions even by using a motion generating machine which can generate a motion 
vector space with two degrees of freedom. 

Fig. 6 shows a coordinate system fixed on the inertia sensor and Fig. 7 shows a 
setup for establishing coincidence between the coordinate axis ZX of the motion generating 
machine and the ZaXa axis of the accelerometer. Then, Fig. 8 shows a setup for 
establishing coincidence between the coordinate axis ZX of the motion generating machine 
and the ZaYa axis of the accelerometer. 

Three sensitivity axes of this accelerometer 1 are denoted as X axis, Y axis and Z 
axis and the Laplace-transformed output signals are denoted by (aox(ja>X *oyQ®>\ &oz (j®))- 
The Laplace-transformed input accelerations are denoted by (aix(ja>), a iy (j(D), aiz (jco)). It is 
not stipulated generally that the input acceleration falls on the vector space (sensitivity 
space) determined by the three sensitivity axes of the accelerometer. When this stipulation 
is accepted, it does not result in supporting a conclusion that the input acceleration is 
regarded as a vector. The accelerometer 1 projects the input acceleration space on the 
sensitivity vector space. 

At this time, the matrix sensitivity is expressed by the matrix of 3 x 3 shown below. 

(Mathematical 22) 

f Cf o o N 

e e e 

Each of Sxpc, Sx, y and expresses the relation of the input signal to the 
accelerometer 1 to the X axis output of the accelerometer 1. Sxx denotes the normal 
sensitivity expressing the relation of the X axis output to the X axis input component, Sx, y 
denotes the cross sensitivity expressing the relation of the X axis output signal to the Y axis 
input component, and Sx, z denotes the cross sensitivity expressing the relation of the X axis 
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output signal to the Z axis input component of the acceleration. Each of S y , x , S y , y and S y ^ 
expresses the relation of the input signal to the Y axis out of the accelerometer 1. Sy* 
denotes the cross sensitivity expressing the relation of the Y axis signal to the X axis input 
component of the accelerometer 1, S y , y denotes the normal sensitivity expressing the 
relation of the Y axis output signal to the Y axis input component of the acceleration, and 
S y , z denotes the cross sensitivity expressing the relation of the Y axis output signal to the Z 
axis input component of the acceleration. 

Each of Sz, x , Sz, y and Sz^ expresses the relation of the input signal to the Z axis 
output of the accelerometer 1. Sz, x denotes the cross sensitivity expressing the relation of 
the Z axis output signal to the X axis input component of the accelerometer 1, Sz, y denotes 
the cross sensitivity expressing the relation of the Z axis output signal to the Y axis input 
component of the accelerometer 1, and Sz, z denotes the normal sensitivity expressing the 
relation of the Z axis output signal to the Z axis input component of the accelerometer 1 . 
At this time, the relation between the input vector and the output vector is expressed by the 
following formula. 

(Mathematical 23) 

S z ^ x {co),S zy {a>),S z%z {a>) 



< a^exp(.;^r) > 



The preceding description is the same as what has been stated in the portion dealing 
with the classification in the third aspect of the invention. Concrete examples of the 
method for defining the sensitivity matrix of a three-axis accelerometer 1 used as an inertia 
sensor for a two-dimensional motion generating machine as a matrix of 3 x 3, namely the 
method for setting up a scene of the calibration, are illustrated in Fig. 6 - Fig. 8. Fig. 6 is a 
perspective view of a coordinate system X, Y, Z fixed on an accelerometer 1 subject to the 
calibration contemplated by this invention, Fig. 7 is a scene set up for establishing 
coincidence between the coordinate axis ZX of the motion generating machine and the 
ZaXa axis of the accelerometer 1 subject to the calibration of this invention, and Fig. 8 is a 
scene set up for establishing coincidence between the coordinate axis ZX of the motion 
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generating machine and the ZaYa axis of the accelerometer 1 subject to the calibration of 
this invention. 

A two-dimensional motion generating machine is only capable of generating two 
independent vector vibration accelerations in a two-dimensional space. So, the vibration 
vector acceleration which can be formed in the setup of Fig. 7 is represented by (a^i, 0, 
aizi)exp(ja>t), (aix2, 0, aiz2)exp(ja)t)and the vibration vector acceleration which can be formed 
in the setup of Fig. 8 is represented by (0, ai y 3, aiz3)exp(jcot). Here, the three vectors (a^i, 
0, aizi), (aix2, 0, a^) and (0, ai y3 , a^) must be independent vectors in a three-dimensional 
space. When these vibration vector accelerations are substituted in Mathematical 23, the 
following three equations, i.e. Mathematical 24, Mathematical 25 and Mathematical 26 are 
established. The terms ao X , ao y and a<>z are complex numbers as the functions of a> which 
takes into account the initial phase of exp(jcot) and the gain. When the three equations, i.e. 
Mathematical 24, Mathematical 25 and Mathematical 26 are unified as simultaneous linear 
equations concerning all the components of the sensitivity matrix, Mathematical 26 is 
derived. 

(Mathematical 24) 



i*^<®>«PC/®*), 



Sy tX (&),S yy (a>) t S yz (a>) 



0 



^r &1 exp( ja>t) 



J 



(Mathematical 25) 

,««2 w«cpt/»oJ [wy^w > 



0 
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(Mathematical 26) 



S yx {a>),S yy {&),S y t {a>) a iy3 expO»0 



(Mathematical 27) 
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< a 0 z3> 



Since the three vectors (a^, 0, a^), (k = 1, 2) and (0, ai y3 , a^) in Mathematical 27 
are linearly independent, it is only plain that the 9 vectors shown in the coefficient 
determinant 

(abdc, 0, aizk , 0, 0, 0, 0, 0, 0) (k = 1, 2), 

(0, 0, 0, aixk , 0, ai2k , 0, 0, 0) (k = 1, 2), 

(0, 0, 0, 0, 0, 0, 0, aixk , 0, aizk) (k = 1, 2), 

(0,a iy3 , ai*, 0 5 0, 0, 0, 0, 0), 

(0, 0, 0, 0, ai y3 , a^, 0, 0, 0) and 

(0, 0, 0, 0, 0, 0, 0, a iy3 , as*) 
are linearly independent. Since the coefficient determinant of Mathematical 27 is not zero, 
it is only logical to conclude that they are capable of being solved without fail. The 
sensitivity Sij (i = x, y, z, j = x, y, z) is a normal sensitivity when the suffixes are equal or a 
cross sensitivity when the suffices are not equal. By scanning the angular frequencies with 
the simultaneous linear equations of Mathematical 27, it is made possible to define the 
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normal sensitivity and the cross sensitivity as the functions of angular frequency. In this 
example, two independent vibration vector accelerations are selected for the setup of Fig. 7 
and one vibration vector acceleration is selected for the setup of Fig. 8. Alternatively, 
one-vibration vector acceleration may be selected for the setup of Fig. 7 and two 
independent vibration vector accelerations may be selected for the setup of Fig. 8. 

Here, the fact that the simultaneous linear equations for deriving the matrix 
sensitivity can be solved will be verified. 

The question whether or not the simultaneous linear equations represented by 
Mathematical 2, Mathematical 8, Mathematical 14, Mathematical 21 and Mathematical 27 
are solvable is very important for the technique disclosed by this invention. Thus, the fact 
that the determinant of N x M-dimension coefficient matrix concerning the sensitivity of an 
inertia sensor having N degrees of freedom on which a motion generating machine having 
M degrees of freedom exerts M independent linear motion vectors is not zero is verified as 
follows. It is assumed that the linearly independent motion vectors imparted by the 
motion generating machine are given by the following formula. 

(Mathematical 28) 
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Since the preceding vectors are linearly independent, the verification of 
Mathematical 29 requires all the coefficients Ai, A 2 , A 3 , Am to be zero. 
(Mathematical 29) 
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What matters resides in the question whether or not the following vector is linearly 
independent under the assumption of Mathematical 29. It is assumed that the linearly 
independent M motion vectors of Mathematical 28 are represented by ai, a 2 , a 3 , a M and 
the zero vector having M O's in a row is represented by 0 having an arrow mark attached 
thereto. The N x M vectors for the investigation regarding the identity of linear 
independence are as follows: 

(Mathematical 30) 
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The fact that all the coefficients are zeroed when the linear combination of the 
vectors shown in Mathematical 30 is set at zero verifies that the vectors are linearly 
independent. Since the renewal of the linear combination merely results in producing N 
formulas conforming with Mathematical 29 and since ai, a 2 , a 3 , ... , a M are linearly 
independent, all the coefficients are proved to be zero and the determinant of coefficients 
matrix evidently is not zeroed. 
(Example 6) 

Example 6 constitutes itself an example of avoiding the installation of an inertia 
sensor 4 for the control of a motion generating machine as illustrated in Fig. 9. Similarly 
to the classification of the first aspect of the invention dealt with in Example 1, the case of 
adopting a one-axis accelerometer as the accelerometer 1 subject to calibration and using a 
three-dimensional motion generating machine will be explained here. 

The operating part of a processor 9 is manipulated to transmit to actuators 5, 6 and 
7 control signals adapted to induce prescribed vibration therein. By applying the 
prescribed signals individually to the actuators 5, 6 and 7 and using one sensitivity axis of 
an accelerometer 1 as X axis, the output signal from the accelerometer 1 is enabled to be 
input to the processor 9. The input which has undergone Laplace transform is represented 
by (ao X (jo>), 0, 0). The letter a> denotes the angular frequencies. The Y axis component 
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and the Z axis component are both zero of course. 

In laser interferometers 2 and 3 for real-time control or measurement, which are 
means to measure the input to the accelerometer 1, the input acceleration to the 
accelerometer 1 is measured through a table 11. The measured input acceleration is input 
to the processor 9. The input which has undergone Laplace transform is represented by 
(aixOoo), aiy(ja>) 5 a^Oo))). It is not stipulated that the input acceleration falls on the 
sensitivity axis of the accelerometer 1. When this stipulation is accepted, it does not result 
in supporting a conclusion that the input acceleration is regarded as a vector. 

At this time, since an inertia sensor 4 for controlling a motion generating machine is 
not installed, no measurement can be made to judge whether or not the actuators 5, 6 and 7 
are generating prescribed motions and consequently the processor 9 cannot emit control 
signals to the actuators 5, 6 and 7. The accelerometer 1 can be calibrated, however, 
because the means 2 and 3 for measuring the input to the accelerometer 1 and the 
processor 9 as a means for measuring the output from the accelerometer 1 and other 
pertinent components are provided. 

Incidentally, the "Laplace transform" mentioned above is intended to ensure 
mathematical strictness of the description of the solution of an equation. It may be 
performed generally by integral transformation or particularly practically by high speed 
Fourier transformation (FFT). It is accomplished so long as the processor 9 is furnished 
with a signal processing function capable of converting the functions in the time domain to 
the functions in the frequency domain. 

Now, the aforementioned processing of matrix sensitivity will be described more 
specifically below. 

The degrees of freedom of a motion in a space total 6, i.e. 3 degrees of freedom of 
a translation^ motion and 3 degrees of freedom of a rotary motion. The largest value of 
the degrees of freedom of an inertia sensor which detects a motion, therefore, is 6 and the 
largest value of the degrees of freedom of a motion generating device is 6. Regarding the 
connection between the degrees of freedom of the inertia sensor subject to calibration and 
the degrees of freedom of the motion which can be generated by the motion generating 
machine, since various combinations can be set, the degrees of freedom will be explained as 
shown in the following table. 
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(Table 1) 





Degree of freedom of motion generation machine (N) 


I 


^ 2 


3 


4 


5 


6 


Degree of 

freedom 

of 

detection 

inertia 

sensor 

(number 

of axis) 

(M) 


1 


This column 
covers 
conventional 
technique . 


The transverse 

sensitivity, 

though 

assumable, 

defies perfect 

definition. 


Explanation 
in (1-3) 


Explanation 
in (1-4) 


Explanation 
in (1-5) 


Explanation 
in (1-6) 


2 


Explained in the case of a biaxial acceleration sensor described herein beiow. 


3 


Explained in the case of a triaxial acceleration sensor described herein below. 


4 


Explained in the case of the general discussion made herein below 


5 


Explained in the case of the general discussion made herein below. 


6 


Explained in the case of the general discussion made herein below 



It is important in any of the cases excepting the case of N = M = 1 that an equation 
can be set in accordance with the following concept. 

(1-3) The case of a one-axis accelerometer (the case of calibrating a one-axis 
accelerometer with a motion generating machine having 3 degrees of freedom) 

The case of handling a one-axis accelerometer as a one-dimensional accelerometer 
and deriving the sensitivity matrix as 1 x 3 matrix will be explained below. 

1) The case of using sinusoidal wave for input acceleration 

It is assumed that the dynamic sensitivity matrix of a one-axis accelerometer at a 
given frequency fi(©i=27rfi) is represented by (Sx^(a>i), Sx, y ((Oi), Sx^(a>i)). The act of 
finding the values of all the components of this matrix constitutes itself the concept of 
calibration. Simultaneous linear equations having these three components as unknown 
quantities are set as shown below. 

(Mathematical 31) 
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Procedure 1 

A vibration is imparted with a given frequency fi(a>i=27rfi) first in the direction of 



-40- 

the vector (a^i, ai y i, aizi). The input acceleration vector signal is represented by (a^i, ai y i, 
aizi)exp(j©t). The output signal of the accelerometer at this time is represented by 
aoxi(a>)exp(ja>t). Since the actual vibration assumes the form of a burst signal and the 
signal of the vibration vector (a^i, ai y i, aizi)exp(ja>t) consequently assumes the form of a 
burst signal, the three components ought to be considered as resulting from extracting the 
fi component through FFT. 

In other words, when the actual vibration signal on the time base is expressed by the 
function of the vector (aixi, ai y i, aizi)Bi(t), the part Bi(t) means the burst signal. It is 
assumed that the result of extracting the fi component by subjecting this burst part to FFT 
is the complex number of Z is i. By the same token, the fi component is extracted through 
FFT because the output signal of the accelerometer at this time is in the form of a burst 
signal. Since the fi component which is obtained at this time is a complex number, the 
division of this component by exp(ja>t) results in ao X i(a>i). 

In other words, the fi component is derived through FFT because the output signal 
of the accelerometer is in the form of burst signal as well. Here, the coefficients 
mentioned above are decided as follows: 

(Mathematical 32) 

Procedure 2 

A vibration is imparted with a given frequency fi(coi=27cfi) first in the direction of 
the vector (a^, ai y2 , a^). The input acceleration vector signal is represented by (a^, ai y2 , 
aiz2)exp(ja>t). The output signal of the accelerometer at this time is represented by 
aox2(co)exp(jcot). Since the actual vibration assumes the form of a burst signal and the 
signal of the vibration vector (aix2, aj y2 , a^expGGOt) consequently assumes the form of a 
burst signal, the three components ought to be considered as resulting from extracting the 
fi component through FFT. By the same token, the fi component is extracted through 
FFT because the output signal of the accelerometer at this time is in the form of a burst 
signal. Since the fi component which is obtained at this time is a complex number, the 
division of this component by exp(ja>t) results in aox2(a>i). 
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When the actual vibration signal on the time base is expressed by the functions of 
the vector (a«2> ai y2 , aiz 2 )B 2 (t), the part B 2 (t) means the burst signal. It is assumed that the 
result of extracting the fi component by subjecting this burst part to FFT is the complex 
number of Z^- By the same token, the fi component is extracted through FFT because 
the output signal of the accelerometer at this time is in the form of a burst signal. Since 
the fi component which is obtained at this time is a complex number, the division of this 
component by exp(ja>t) results in a^coi). In other words, the fi component is derived 
through FFT because the output signal of the accelerometer is in the form of burst signal as 
well. It is assumed that this signal is a complex number of Z os2 . Here, the coefficients 
mentioned above are decided as follows: 

(Mathematical 33) 

C 2,1 = a ix2 C Z2 = a iy2 C 2,3 = a iz2 d 2 = Q ox2 = Z os2 lZ is2 
Procedure 3 

A vibration is imparted with a given frequency fi((Oi=27tfi) first in the direction of 
the vector (a^, ai y3 , a^). The input acceleration vector signal is represented by (a^, ai y3 , 
ai23)exp(j(Dt). The output signal of the accelerometer at this time is represented by 
a ox3 (co)exp(j(Dt). Since the actual vibration assumes the form of a burst signal and the 
signal of the vibration vector (a^, aj y3 , aiz3)exp(ja>t) consequently assumes the form of a 
burst signal, the three components ought to be considered as resulting from extracting the 
fi component through FFT. By the same token, the fi component is extracted through 
FFT because the output signal of the accelerometer at this time is in the form of a burst 
signal. Since the fi component which is obtained at this time is a complex number, the 
division of this component by exp(jcot) results in a^cDi). 

When the actual vibration signal on the time base is expressed by the function of the 
vector (aix3, aj y3 , a^B^t), the part B 3 (t) means the burst signal. It is assumed that the 
result of extracting the fi component by subjecting this burst part to FFT is the complex 
number of Zj s3 . By the same token, the fi component is extracted through FFT because 
the output signal of the accelerometer at this time is in the form of a burst signal. Since 
the fi component which is obtained at this time is a complex number, the division of this 
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component by exp(ja>t) results in aox3(a>i). In other words, the fi component is derived 
through FFT because the output signal of the accelerometer is in the form of burst signal as 
well. This signal is a complex number of Z 0S 3. Here, the coefficients mentioned above 
are decided as follows: 

(Mathematical 34) 

C 3,1 = *«3 C 3,2 = a iy3 C 3.3 = a iz3 d 3 = a ax3 = Z os3 f Z h3 
Procedure 4 

Consequently, the coefficient matrix of the aforementioned simultaneous linear 
equations has been decided and the vectors on the right part have been decided as well. 
The aforementioned simultaneous linear equations are solved without fail because the 
vectors in the three directions are selected in linear independence. Subsequently, the 
matrix sensitivity in a different frequency can be found by setting the frequency of vibration 
at another value and carrying out the same calculation. 

2) The case of using a random waveform for input acceleration 

The direction vector of random vibration is selected in the direction of linear 
independence. This direction is represented as shown below: 

(aixi, ai y i, aizi), (a^, ai y2 , a^) and (a^, ai y3 , a^). 
The frequency band of the random vibration signal in these directions is represented as 
follows: 

[fmin fmax]- 

It is assumed that the dynamic sensitivity matrix of the one-axis accelerometer at a given 
frequency fi(a>i=27tfi) is represented by (Sx^(odi), S^cdi), Sx, z ((Oi)). The act of finding the 
values of all the components of this matrix constitutes itself the concept of correction. 
Simultaneous linear equations having these three components as unknown quantities are set 
as shown below. 
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(Mathematical 35) 
c 2,1 c 2,2 c 2,3 









= d 2 







Procedure 1 

A random vibration is imparted in the direction of (a«i, ai y i, aizi). The output 
signal of the accelerometer which is obtained at this time is represented by the function of 
time f<,i(t). The input signal vector is represented by (aixi, ai y i, aizi)fii(t). The random 
function is fu(t). Here, f u (t) and £>i(t) are subject to FFT. As a result, the data shown in 
the following table is obtained. 



(Table 2) 



Frequency 


Result of FFT of random input signal 


Result of FFT of output signal of 
acceleration sensor 
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max 


mat 



Here, the aforementioned coefficients are decided as follows: 



(Mathematical 36) 

C U =a ixl C M= a iyi C ^= a izl <ti= a <Ki = Z o,i/ Z iA 



Procedure 2 

A random vibration is imparted in the direction of (a^, ai y2 , a^). The output 
signal of the accelerometer which is obtained at this time is represented by the function of 
time £,2(t). The input signal vector is represented by (a^, aj y2 , a^foO). The random 
function is f i2 (t). Here, f i2 (t) and f o2 (t) are subject to FFT. As a result, the data shown in 
the following table is obtained. 
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(Table 3) 



Frequency 


Remit of FFT of random input signal 


Result of FFT of output signal of 
acceleration sensor 
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www 
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Here, the aforementioned coefficients are decided as follows: 
(Mathematical 37) 

^2,1 = fl ix2 ^2,2 ~ a iy2 ^2,3 = Q iz2 ^2 ~ a ox2 ~ % o,2 /%i 9 2 
Procedure 3 

A random vibration imparted in the direction of (a«3, ai y3> a^). The output signal 
of the accelerometer which is obtained at this time is represented by the function of time 
f 0 3(t). The input signal vector is represented by (a^, ai y3 , a^feCt). The random function 
is fi 3 (t). Here, fl 3 (t) and f G 3(t) are subject to FFT. As a result, the data shown in the 
following table is obtained. 



(Table 4) 



Frequency 


Result of FFT of random input signal 


Result of FFT of output signal of 
acceleration sensor 


4inOmin) 


min 










<*>! 












^max( ^zrux) 


max 





Here, the aforementioned coefficients are decided as follows: 
(Mathematical 38) 

^3,1 = a ix3 C\2 — a iyl Q,3 = a iz3 ^3 = a ox3 = ^o 3 3 / ^ i,3 

Procedure 4 

Consequently, the coefficient matrix of the aforementioned simultaneous linear 
equations has been decided and the vectors on the right part have been decided as well. 
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The aforementioned simultaneous linear equations are solved without fail because the 
vectors in the three directions are selected in linear independence. 

Subsequently, the matrix sensitivity in a different frequency can be found by 
selecting another frequency in the table obtained by spectral analysis and performing the 
same calculation. By performing this calculation on all the frequencies obtained by 
spectral analysis, it is made possible to find the matrix sensitivity at [f„un fmax] by the 
frequency resolving power which is decided by the setting of FFT to be used in that case. 

3) The case of using pulse for input acceleration 

The direction vector of input acceleration with pulses is selected in the direction of 
linear independence. This direction is represented as shown below: 

(aixi, ai y i, a^i), (a*2, ai y2 , a^) and (a^, ai y3 , a^). 
The frequency band of the pulse input acceleration signal in these directions is represented 
as follows: 

[fmin dux]- 

It is assumed that the dynamic sensitivity matrix of the one-axis accelerometer at a given 
frequency fi(a>r=27cfi) is represented by (Sx, x (o3i), ^Sx, y (a>i), Sx^(a>i)). The act of finding 
the values of all the components of this matrix constitutes itself the concept of calibration 
contemplated by this invention. Simultaneous linear equations having these three 
components as unknown quantities are set as shown below. 
(Mathematical 39) 



c 2,1 C Z2 c 2,3 









= d 2 







Procedure 1 

A pulse input acceleration is imparted in the direction of (aixi, ai y i, aizi). The 
output signal of the accelerometer which is obtained at this time is represented by the 
function of time f 0 i(t). The input signal vector is represented by (aixi, aj y i, aizi)fii(t). The 
pulse function is fh(t). Here, f u (t) and f G i(t) are subject to FFT. As a result, the data 
shown in the following table is obtained. 
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(Table 5) 



Frequency 


Result of FFT of pulse acceleration 
input signal 


Result of FFT of output signal of 
acceleration sensor 








♦ 




* 


CO, 
















aast 



Here, the aforementioned coefficients are decided as follows: 
(Mathematical 40) 

C l,l =a ix\ C i 9 2 =a iy\ C \J= a iz\ d \= a ox\= Z oA/ Z u\ 
Procedure 2 

A random vibration is imparted in the direction of (aix2, ai y2 , a^). The output 
signal of the accelerometer which is obtained at this time is represented by the function of 
time f 0 2(t). The input signal vector is represented by (a^a, ai y2 , a^faCt). The random 
function is f i2 (t). Here, f i2 (t) and f o2 (t) are subject to FFT. As a result, the data shown in 
the following table is obtained. 

(Table 6) 



Frequency 


Result of FFT of pulse acceleration 
input signal 


Result of FFT of output signal of 
acceleration sensor 


^nan(^inin) 


^i. rain. 


■^ta. min 



























Here, the aforementioned coefficients are decided as follows: 
(Mathematical 41) 

C 2,l = a ix2 C 2a= a iy2 C 2,3=^2 d 2 = <* o*2 = Z o,2 / Z 1,2 
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Procedure 3 

A random vibration is imparted in the direction of (aix3, ai y3 , a^). The output 
signal of the accelerometer which is obtained at this time is represented by the function of 
time f 0 3(t). The input signal vector is represented by (aix3, ai y3 , a^feCt). The random 
function is fi3(t). Here, f i3 (t) and f o3 (t) are subject to FFT. As a result, the data shown in 
the following table is obtained. 

(Table 7) 



Frequency 


Result of FFT of pulse acceleration 
input signal 


Result of FFT of output signal of 
; acceleration sensor 
















2^ 


Zo.3 















Here, the aforementioned coefficients are decided as follows: 
(Mathematical 42) 

^3,1 =C *ix3 =a iy3 ^3,3 = G iz3 ^3 =Q ox3 = % o,3 / ^ i y 3 

Procedure 4 

Consequently, the coefficient matrix of the aforementioned simultaneous linear 
equations has been decided and the vectors on the right part have been decided as well. 
The aforementioned simultaneous linear equations are solved without fail because the 
vectors in the three directions are selected in linear independence. Subsequently, the 
matrix sensitivity in a different frequency can be found by selecting another frequency in the 
table obtained by spectral analysis and performing the same calculation. By performing 
this calculation on all the frequencies obtained by spectral analysis, it is made possible to 
find the matrix sensitivity at [fmin fmax] by the frequency resolving power which is decided by 
the setting of FFT to be used in that case. 

The method for finding the matrix sensitivity of the one-axis accelerometer by using 
a three-axis motion generating machine has been described. Actually, the case of finding 
the matrix sensitivity of a one-axis accelerometer by using a four-axis motion generating 
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machine, the case of finding the matrix sensitivity of a one-axis accelerometer by using a 
five-axis motion generating machine, and the case of finding the matrix sensitivity of a 
one-axis accelerometer by using a six-axis motion generating machine are conceivable. 
These cases may be severally coped with by using the following approaches. 

(1-4) The case of finding matrix sensitivity of one-axis accelerometer by using a 
four-axis motion generating machine 

(Note: No particular need is found for designating which of the four axes are used 
for translational acceleration or for rotational acceleration.) The case of handing a 
one-axis accelerometer as a one-dimensional accelerometer deriving the sensitivity matrix 
as a matrix of 1 x 4 will be explained below. 

1) The case of using sinusoidal wave for input acceleration 

The output axis of an inertia sensor is taken as the X axis. The axes of motion 
generated by ef a motion generating machine are taken as X, Y, Z and a. The dynamic 
sensitivity matrix of a one-axis accelerometer at a given frequency fi(a>i=27ifi) is 
represented by (Sx^((Di) ? Sx,y(a>i), Sx^(coi), S^an)). The act of finding all the components 
of this matrix constitutes itself the concept of calibration. Simultaneous linear equations 
using these four variables as unknown quantities are set as follows. 

(Mathematical 43) 





C|_2 




<v 






f*0 








C 2.4 








C 3k1 














, C 4.1 


C42 


C AZ 











Since the motions which can be generated by a motion generating machine can form 
a four-dimensional space, it is important that four independent vectors be selected. 
Procedure 1 

Since the motions which can be generated by the motion generating machine are 
capable of forming a four-dimensional space and allowing selection of four independent 
vectors, Procedure 1 first selects the first vector and imparts a sinusoidal wave burst signal 
having a frequency fi(a>i=27rfi) in that direction. The vector selected from among the 
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independent vectors is represented by (a^i, &i y u &izi> aiai). The input vector acceleration at 
this time is represented by (a^i, ai y i, aizi, a io a)Bi(t). At this time, the output signal of the 
inertia sensor assumes the form of a sinusoidal wave burst signal. It is assumed that the 
spectral analysis finds the component of the frequency fi(cc>i=27tfi) to be Z 0 , s m. By this 
procedure, the first row of the coefficient matrix of the simultaneous linear equations 
(Mathematical 43) and the first row of the constant vector are found as shown below. 
(Mathematical 44) 




Procedure 2 

Since the motions which can be generated by the motion generating machine are 
capable of forming a four-dimensional space and allowing selection of four independent 
vectors, Procedure 2 first selects the second vector and imparts a sinusoidal wave burst 
signal having a frequency fi(a>i=27cfi) in that direction. The vector selected from among 
the independent vectors is represented by (a^, ai y2 , a^, a^). The input vector 
acceleration at this time is represented by (a^, a iy2 , a^, aia2)B 2 (t). At this time, the output 
signal of the inertia sensor assumes the form of a sinusoidal wave burst signal. It is 
assumed that the spectral analysis finds the component of the frequency fi(a>i=27ifi) to be 
Zi,sb2- By this procedure, the second row of the coefficient matrix of the simultaneous 
linear equations (Mathematical 43) and the second row of the constant vector are found as 
shown below. 

(Mathematical 45) 

c 2,1 = *fr2 C 2 % 2 = °iy2 C 23 = ^z2 C ZA = ^2 




Procedure 3 

Since the motions which can be generated by the motion generating machine are 
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capable of forming a four-dimensional space and allowing selection of four independent 
vectors, Procedure 3 first selects the third vector and imparts a sinusoidal wave burst signal 
having a frequency fi(coi=27cfi) in that direction. The vector selected from among the 
independent vectors is represented by (a^, ai y3 , a^, ai^). The input vector acceleration at 
this time is represented by (a^, ai y3 , a^, a ia3 )B 3 (t). At this time, the output signal of the 
inertia sensor assumes the form of a sinusoidal wave burst signal. It is assumed that the 
spectral analysis finds the component of the frequency fi(a>i=27cfi) to be Zi, sb3 . By this 
procedure, the third row of the coefficient matrix of the simultaneous linear equations 
(Mathematical 43) and the third row of the constant vector are found as shown below. 
(Mathematical 46) 




Procedure 4 

Since the motions which can be generated by the motion generating machine are 
capable of forming a four-dimensional space and allowing selection of four independent 
vectors, Procedure 4 4- first selects the fourth vector and imparts a sinusoidal wave burst 
signal having a frequency fi(a>i=27tfi) in that direction. The vector selected from among 
the independent vectors is represented by (a^, ai y4 , a^, aic*). The input vector 
acceleration at this time is represented by (a^, ai y4 , a^, aio4)B 4 (t). At this time, the output 
signal of the inertia sensor assumes the form of a sinusoidal wave burst signal. It is 
assumed that the spectral analysis finds the component of the frequency fi(a>i=27cfi) to be 
Z 0 ,sb4- By this procedure, the fourth row of the coefficient matrix of the simultaneous 
linear equations (Mathematical 43) and the fourth row of the constant vector are found as 
shown below. 
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(Mathematical 47) 

C AA =a ix4 C 4.2 = a iyA C AZ =a izA C 4A =a toA 




(1-5) The case of finding the matrix sensitivity of one-axis accelerometer by using 
five-axis motion generating machine 

(Note: No particular need is found for designating which of the five axes are used 
for translational acceleration or for rotational acceleration.) The case of handing a 
one-axis accelerometer as a one-dimensional accelerometer deriving the sensitivity matrix 
as a matrix of 1 x 5 will be explained below. 

1) The case of using a sinusoidal wave for input acceleration 

The output axis of an inertia sensor is taken as the X axis. The axes of motion 
generated by a motion generating machine are taken as X, Y, Z, a, p. The dynamic 
sensitivity matrix of a one-axis accelerometer at a given frequency fi(a>i=27tfi) is 
represented by (Sx, x (coi), S^coi), Sx, 2 (a>i), Sx, a (coi), Sxp(a>i)). The act of finding all the 
components of this matrix constitutes itself the concept of calibration. Simultaneous 
linear equations using these five variables as unknown quantities are set as follows: 



(Mathematical 48) 
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(1-6) The case of finding matrix sensitivity of one- axis accelerometer by using 
six-axis motion generating machine (Note: No particular need is found for designating 
which of the five axes are used for translational acceleration or for rotational acceleration.) 

The case of handing a one-axis accelerometer as a one-dimensional accelerometer 
deriving the sensitivity matrix as a matrix of 1 x 6 will be explained below. 
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1) The case of using a sinusoidal wave for input acceleration 

The output axis of an inertia sensor is taken as the X axis. The axes of motion 
generated by a motion generating machine are taken as X, Y, Z, a, P, y. The dynamic 
sensitivity matrix of a one-axis accelerometer at a given frequency fi(a>i=27ifi) is 
represented by (S^coi), Sx,y(coi), Sx, z (g)i), Sx, a (a>i), Sx,p(a>i)}, Sx, Y (a>i)). The act of finding 
all the components of this matrix constitutes itself the concept of calibration. 
Simultaneous linear equations using these six variables as unknown quantities are set as 
follows. 

(Mathematical 49) 

* 
dA 

<** 

The case of two-axis accelerometer 

The sensitive axes of a two-axis accelerometer are taken as X axis and Y axis. 
The Laplace-transformed output signal is taken as (ao X (a>), aoxy(co), 0). The letter a> 
denotes the angular frequencies. The Z component is zero because the accelerometer is a 
two-axis sensor. The Laplace-transformed input acceleration is represented by (aix(co), 
aiy(a>), aiz(co)). It is not stipulated that the input acceleration falls on the plane (sensitivity 
plane) which is decided by two sensitivity axes of the accelerometer. When this 
stipulation is accepted, it does not result in supporting a conclusion that the input 
acceleration is regarded as a vector. At this time, the matrix sensitivity is expressed by a 
matrix of 2 x 3 shown below. 
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( 



(Mathematical 50) 
O C C \ 

e e o 
a y 9 x^y 9 y^y 9 zJ 



Incidentally, the act of finding all the components of this sensitivity matrix 
constitutes itself the concept of calibration contemplated by this invention. Sx, x , Sx, y , S^, 
Sy^c, Sy, y and S y , z are found from the following simultaneous linear equations. It is 
provided, however, that the coefficient matrix of the simultaneous linear equations and the 
vector on the right part are more or less varied when the vibration is made with a sinusoidal 
wave, when the vibration is made with a random waveform, and when the vibration is made 
with a pulse waveform. 

(Mathematical 51) 







C \3 


C t4 


^15 












<%2 


C23 
















0*2 


C33 


3m 










<*3 




C42 






C4.5 








dA 




°52 








^6 






** 




Qt2 


^63 


^4 


<%6 






<*6, 



The procedures for deriving the coefficient Qj (i, j = 1, 6), dj (i = 1, 6) from 
the preceding formula will be shown below. 

1) The case of using a sinusoidal wave for vibration 

When the vibration is made with a sinusoidal wave, theoretically the relation 
between the input acceleration vector (a^ a iy , aiz)exp(ja>t) and the output acceleration 
vector (aox, ao y )exp(ja>t) is as shown below. 
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(Mathematical 52) 



^exp( jmt) j 



Procedure 1 

Theoretically, a vibration is imparted first in the direction (a^, ajyi, a^) of the vector 
at a given frequency fi(coi=27rfi). The input acceleration vector signal is expressed by (a^i, 
aiyi, aizi)exp(jo>it). The X axis output signal and the Y axis output signal of the 
accelerometer are respectively represented by ao X i(a>i)exp(ja>it) and ao y i(a>i)exp(ja>it). 
Since the actual vibration is made with a burst signal and consequently the signal (a^i, ai y i, 
aizi)exp(jcL)t) of the vibration vector is also a burst signal, the three components ought to be 
considered as resulting from the extraction of the fi component through FFT. In other 
words, when the actual vibration signal is expressed on the time base by the function of the 
vector (aixi, ai y i, a^Bi (t), the part Bi(t) means the burst signal. It is assumed that the 
complex number Z^i results from extracting the fi component from the burst part through 



Since the output signals (B ox i(t), B oy i(t)) of the accelerometer in this case are also 
burst signals, the fi component is extracted likewise through FFT. The fi component 
which is obtained at this time is a complex number, the division of this component by Z^i 
results in (aoxi(<Oi), ao y i(a>i)). As a result, the following formula is established for the 
frequency component 1 . 

(Mathematical 53) 



When the coefficient of Mathematical 5 1 is derived by unifying this formula into 
equations concerning (Sx, x , S^y, Sx^, S y , x , S y , y , S y , z ), the following formula is established. 



FFT. 
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(Mathematical 54) 



<*1=*a*1<»l> * = Vl^) 



Procedure 2 

Theoretically, a vibration is imparted first in the direction ^(aba, ai y 2, a^) of the 
vector at a given frequency fi(a>i=27tfi). The input acceleration vector signal is expressed 
by (a^, aj y2 , ai22)exp(jcDit). The X axis output signal and the Y axis output signal of the 
accelerometer are respectively represented by aox2(©i)exp(jcoit) and ao y 2(o)i)exp(ja)it). 

Since the actual vibration is made with a burst signal and consequently the signal 
(aM, ai y2 , aiz2)exp(ja>t) of the theoretical vibration vector is also a burst signal, the three 
components ought to be considered as resulting from the extraction of the fi component 
through FFT. In other words, when the actual vibration signal is expressed on the time 
base by the function of the vector (a^, a^, ai^Bift), the part B 2 (t) means the burst signal. 
It is assumed that the complex number results from extracting the fi component from 
the burst part through FFT. 

Since the output signals (B OX 2(t), B oy2 (t)) of the accelerometer in this case are also 
burst signals, the fi component is extracted likewise through FFT. The fi component 
which is obtained at this time is a complex number, the division of this component by 
results in (aox2(a>i), ^((fli)). As a result, the following formula is established in 
Procedure 2 for the frequency component 1 . 

(Mathematical 55) 
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When the coefficient of Mathematical 5 1 is derived by unifying this formula into 
equations concerning (Sx, x , Sx, y , Sx, Zj S y ^, S y , y , S y , z ), the following formula is established. 
(Mathematical 56) 

C&.1=«*2 c i2=%2 C« = %2 C a4 =0 03.5 = 0 C^ e = 0 
Procedure 3 

Theoretically, a vibration is imparted first in the direction (a^, ai y3 , a^) of the 
vector at a given frequency fi(a>i=27tfi). The input acceleration vector signal is expressed 
by (a^, ai y3 , ai z3 )exp(j(Oit). The X axis output signal and the Y axis output signal of the 
accelerometer are respectively represented by aox3(coi)exp(ja)it) and ao y3 (a>i)exp(jcoit). 
Since the actual vibration is made with a burst signal and consequently the signal (a^, aj y3 , 
aiz3)exp(ja>it) of the theoretical vibration vector is also a burst signal, the three components 
ought to be considered as resulting from the extraction of the fi component through FFT. 
In other words, when the actual vibration signal is expressed on the time base by the 
function of the vector (a**, a iy3 , a iz3 )B 3 (t), the part B 3 (t) means the burst signal. It is 
assumed that the complex number Z is3 results from extracting the fi component from the 
burst part through FFT. 

Since the output signals (B ox3 (t), B oy3 (t)) of the accelerometer in this case are also 
burst signals, the fi component is extracted likewise through FFT. The ft component 
which is obtained at this time is a complex number, the division of this component by Zi s3 
results in (aox3(o>i), ao y3 (coi)). As a result, the following formula is established in 
Procedure 3 for the frequency component 1 . 

(Mathematical 57) 

C 3J = a t* C 5.2 = a *3 C 3J> C 3A =° C 5.5 =° C 5.6 =° 

C w =0 C w =0 C w =0 C M =a w C 6w5 =a w C M =a M 

d 5 =fl«3(^l) d S = 0 oy3(^l) 
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When the coefficient of Mathematical 5 1 is derived by unifying this formula into 
equations concerning (S*,*, Sx, y , S^, S y , x , S y , y , S y ^), the following formula is established. 
(Mathematical 58) 
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The coeflScient matrix of Mathematical 51 and the coefficient of the right part have 
been decided till Procedure 3. The following formula results from rewriting the 
coefficient matrix. 



(Mathematical 59) 
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When the input vectors (a^i, a iy i, a^i), (a^, a iy2 , a^), (a^, ai y3 , a^) in 
Mathematical 59 are selected in linear independence, the coefficient determinant is not zero. 
Thus, the simultaneous linear equations are solved without fail. The reason can be 
explained by the following transformation. 
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(Mathematical 60) 
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The exchange of the second row and the third row in this formula, results in the 
following formula. 

(Mathematical 61) 
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Here, the exchange of the third row and the fifth row in this formula, results in the 
following formula. 

(Mathematical 62) 
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Further, the exchange of the fourth row and the fifth row in this formula, results in 
the following formula. 

(Mathematical 63) 
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Procedure 4 

Consequently, the coefficient matrix of the aforementioned simultaneous linear 
equations has been decided and the vectors on the right part have been decided as well 
The aforementioned simultaneous linear equations are solved without fail because the 
vectors in the three directions are selected in linear independence. Subsequently, the 
matrix sensitivity in a different frequency can be found by selecting another frequency in the 
table obtained by spectral analysis and performing the same calculation. Specifically, the 
frequency is scanned in a necessary frequency band. 

2) The case of using a random waveform for input acceleration 

The direction of random vibration is selected in a direction of linear independence. 
This direction is represented by (a^i, ai y i, aizi), (a^, ai y2 , a^)? and (aix3, aj y3 , a^). It is 
assumed that the frequency band of random vibration signals in these directions is 
expressed by [fmin f^I The sensitivity axes of a two-axis accelerometer are denoted by X 
axis and Y axis. The Laplace-transformed output signals are represented by (aox(co), 
aoy(to), 0). The Z component is zero because the accelerometer is a two-axis sensor. 
The Laplace-transformed input acceleration is represented by (a^co), a iy (a>), a^©)). 

It is not stipulated that the input acceleration falls on the plane (sensitivity plane) 
which is decided by two sensitivity axes of the accelerometer. When this stipulation is 
accepted, it does not result in supporting a conclusion that the input acceleration is 
regarded as a vector. At this time, the matrix sensitivity is expressed by a matrix of 2 x 3 
shown below. 
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It is assumed that the dynamic sensitivity matrix of the two-axis accelerometer in a 
given frequency fi(a>i=27tfi) is expressed by the following matrix. 
(Mathematical 64) 

^x t x 9 ^x M y^x 9 z 

The act of finding all the components of this matrix constitutes itself the concept of 
calibration. The simultaneous linear equations concerning these six unknown quantities 
are set as follows: 

(Mathematical 65) 
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The coefficient matrix and the constant term of the right part in Mathematical 65 in 
the case of random vibration are found by the following procedure. 
Procedure 1 

A random vibration is imparted in the (aixi, ajyi, aizi) direction. The X axis output 
signal and the Y axis output signal of the accelerometer obtained at this time are 
represented by the functions of time, (aoix(t), a<>iy(t)). The input signal vectors are 
represented by (a^i, ajyi, a^Of^t). The random function is fu(t). Here, fu(t), f 0 i x (t), and 
foiy(t) are subject to FFT. As a result, the data shown in the following table is obtained 
regarding the individual frequency components. 
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(Table 8) 



Frequency 


Result of FFT of random 
acceleration input signal 


Result of FFT of output 
signal iux(t) of acceleration 
sensor 


Result of FFT of output 
signal ioiy(t) of acceleration 
sensor 






■^b, mm. x 


^b, itrin. v 












Zli.. 


Z*. 1. Lx 
















Zb.m«.x 


^3. HUM. y i 



From this table, part of the coefficient of Mathematical 65 is found as follows: 
(Mathematical 66) 

C U = *«1 C l,2 = a iy\ C l,3 = «*1 C l,4 = 0 C l,5 = 0 C l,6 = 0 
C 2,l =0 C 2,2 =0 C 2,3 =0 C 2,4 = «*1 C 2,5 = "frl C 2,6 = 

/ ^i,U / ^i,U 

Procedure 2 

A random vibration is imparted in the (a^, ai y2 , a^) direction. The X axis output 
signal and the Y axis output signal of the accelerometer obtained at this time are 
represented by the functions of time, (a^t), ao2y(t)). The input signal vectors are 
represented by (a^, ai y2 , a^faCt). The random function is f i2 (t). Here, f i2 (t), £>2x(t), f 0 2 y (t) 
are subject to FFT. As a result, the data shown in the following table is obtained 
regarding the individual frequency components. 



(Table 9) 



Frequency 


Result of FFT of random 
acceleration input signal 

m 


Result of FFT of output 
si^tal 42x0) of acceleration 
sensor 


Result of FFT of output 
signal j$2y(t) of acceleration 
sensor 
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atria, x 
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From this table, part of the coefficient of Mathematical 65 is found as follows: 
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(Mathematical 67) 

Q,l = a ix2 C 3,2 = <Ziy2 Q,3 = <*iz1 Q,4 = 0 Q,5 = 0 Q,6 = <> 
C 4,l =0 C 4,2 =0 C 4,3 =0 C 4,4 =**2 C 4,5 = ^>2 C 4,6 = * M 

Procedure 3 

A random vibration is imparted in the (aix3, ai y3 , a^) direction. The X axis output 
signal and the Y axis output signal of the accelerometer obtained at this time are 
represented by the functions of time, (ao3x(t), ao3y(t)). The input signal vectors are 
represented by (aix3, aj y3 , a i23 )fi3(t). The random function is feO). Here, fi 3 (t), £>3x(t), £>3y(t) 
are subject to FFT. As a result, the data shown in the following table is obtained 
regarding the individual frequency components. 

(Table 10) 



Frequency 


Remit of FFT of random 
acceleration input signal 


Result of FFT of output 
signal i^rlt) of acceleration 
sensor 


Result of FFT of output 
signal 4*y(t) of acceleration 
sensor 




































2l M 







From this table, part of the coefficient of Mathematical 65 is found as follows: 
(Mathematical 68) 

Q,l = <*ix3 C 5 t 2 = a iy3 ^5,3 = a iz3 C 5 9 4 =0 Q,5 =0 Q,6 = 0 
C 6,l = 0 C 6,2 = 0 C 6,3 = 0 C 6,4 = **3 C 6,5 = C 6,6 = */z3 

/ Z /,l,3 / Z /,U3 

Procedure 4 

Consequently, the coefficient matrix of the aforementioned simultaneous linear 
equations (Mathematical 65) has been decided and the vectors on the right part have been 
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decided as well. The aforementioned simultaneous linear equations are solved without fail 
because the vectors in the three directions are selected in linear independence. 
Subsequently, the matrix sensitivity in a different frequency can be found by setting the 
vibration frequency at a different value and performing the same calculation. Specifically, 
by scanning the frequency in a necessary frequency band [fmin fmax] and forming equations 
with the individual frequency, it is made possible to find the matrix sensitivity with the 
individual frequency. 

3) The case of using a pulse waveform for input acceleration 

The direction vector of pulse input acceleration is selected in the direction of linear 
independence. This direction is represented as follows: 

(aixi, a iy i, a^i), (a^, a iy2 , a^) and (a^, a iy3 , a^). 
It is assumed that the frequency band of the pulse input acceleration in these directions is 
expressed as follows: 

[fmin fmax]- 

Then, the dynamic sensitivity matrix of the one-axis accelerometer in a given 
frequency fi(a>i=27cfi) is expressed as follows: 
(Mathematical 69) 



^x,x^x t y^x f z 
<? e c 

V y*x* y*y 1 y%z J 



The act of finding all the components of this matrix constitutes itself the concept of 
calibration contemplated by this invention. The simultaneous linear equations having 
these six unknown quantities are set as follows: 
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(Mathematical 70) 
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When the input acceleration is made with a pulse waveform, the coefficients Cy (i = 
1 — 6, j = 1 ••■ 6), di (i = 1—6) of the aforementioned simultaneous linear equations are 
decided by the following procedures. 
Procedure 1 

A pulse acceleration is imparted in the (a^i, ai y i, a^i) direction. The X axis output 
signal and the Y axis output signal of the accelerometer obtained at this time are 
represented by the functions of time, foix(t) and Ciy(t). The input signal vectors are 
represented by (a^i, a iy i, a^Of^t). The pulse function is fn(t). Here, fn(t), £>i x (t), £>i y (t) 
are subject to FFT. As a result, the data shown in the following table is obtained. 



(Table 11) 



Frequency 


Result of FFT of pulse 
acceleration input signal 


Result of FFI of output 
signal of acceleration 
sensor 


Resuh of FFI of output 
signal Aiy(t) of acceleration 
sensor 


"^tun( ^ mm) 






rrria.v 














Zb.LL. 
















■^o. max. x 





By using the results of this table, it is made possible to find part of the coefficients 
of the simultaneous linear equations as follows: 
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(Mathematical 71) 

Q.1 = Q ix\ C l,2 = "ryl C l,3 = *<rl Q.4 = 0 C l,5 = 0 Q,6 = 0 
C 2,l = 0 C 2,2 = 0 C 2,3 = 0 C 2,4 = C 2,5 = %\ C 2,6 = #1 




Procedure 2 

A pulse acceleration is imparted in the (a^, ai y2 , a^) direction. The X axis output 
signal and the Y axis output signal of the accelerometer obtained at this time are 
represented by the functions of time, f C 2x(t) and f<,2y(t). The input signal vectors are 
represented by (a^, a iy 2, a^O^t). The pulse function is f i2 (t). Here, f i2 (t), f o2x (t), £>2y(t) 
are subject to FFT. As a result, the data shown in the following table is obtained. 

(Table 12) 



Frequency 


Result of FFT of pulse 
acceleration input signal 


Result of FFT of output 
signal ^x(t) of acceleration 
sensor 


Result of FFT of output 
signal JctyAt) of acceleration 
sensor 






^™ 














A 1.2 























By using the results of this table, it is made possible to find part of the coefficients 
of the simultaneous linear equations as follows: 
(Mathematical 72) 

C 3,l = a ix2 C 3,2 = «*2 C 3,3 = «»2 Q,4 = 0 Q,5 = 0 Q,6 = 0 
C 4,l =0 C 4,2 =0 C 4,3 =0 C 4.4 = °*2 C 4.5 = ^2 C 4.6 = fl *2 

d > - A,,, d ' ~ A,, 2 

Procedure 3 

A pulse acceleration is imparted in the (a^, aj y3 , a^) direction. The X axis output 
signal and the Y axis output signal of the accelerometer obtained at this time are 
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represented by the functions of time, f 0 3x(t) and £> 3y (t). The input signal vectors are 
represented by (a^, ai y3 , aiz3)fi3(t). The pulse function is f i3 (t). Here, fi3(t), f C 3x(t), f 0 3 y (t) 
are subject to FFT. As a result, the data shown in the following table is obtained. 



(Table 13) 



Frequency 


Result of FFT of pulse 
acceleration input signal 


Result of FFT of output 
signal ifl3 X (t) of acceleration 
sensor 


Result of FFT of output 
signal ^»v(t) of acceleration 
sensor 


^rninC ^ nun) 


min 


min. x 












to, 


^S, 1.3 


Za. 1.3.x 





















By using the results of this table, it is made possible to find part of the coefficients 
of the simultaneous linear equations (Mathematical 70) as follows: 
(Mathematical 73) 

C 5.1 = a ixl C S.2= a iy3 C 5.3 = <*izi C SA =° C 5,S =° C 5.6 = 0 
C 6.1 =0 C 6,2 =0 C 6.3 =0 C 6,4 = ^1x3 C 6.5 =^3 Q.6 =^3 

/ ^/.U 

Procedure 4 

Consequently, the coefficient matrix of the aforementioned simultaneous linear 
equations (Mathematical 70) has been decided and the vectors on the right part have been 
decided as well. The aforementioned simultaneous linear equations are solved without fail 
because the vectors in the three directions are selected in linear independence. 
Subsequently, the matrix sensitivity in a different frequency can be found by setting the 
vibration frequency at a different value and performing the same calculation. Specifically, 
by scanning the frequency in a necessary frequency band [fmin fmax] and forming equations 
with the individual frequencies, it is made possible to find the matrix sensitivity with the 
individual frequencies. 

The case of three-axis accelerometers 

The sensitivity axes of a three-axis accelerometer are represented respectively by X 
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axis, Y axis and Z axis. The Laplace-transformed output signals are represented by 
(aox(co), aoy(a>), ao z (co)). The letter co denotes the angular frequency. The 
Laplace-transformed input acceleration is represented by (aix(a)), ai y (co), a^oo)). At this 
time, the matrix sensitivity is expressed by a matrix of 3 x 3 shown below. The relation 
between the output signal vector and the input vector signal is expressed by Mathematical 
74. 



(Mathematical 74) 





f 








V 



'xje 



'aw 



V / ^ 



S 



SJ5 A 



Incidentally, the act of finding all the components of this sensitivity matrix 
constitutes itself the concept of calibration contemplated by this invention. Sx, x , Sx, y , Sx^, 
Sy^c, Sy, y , S y>z , Sz^, S^y and are found from the following simultaneous linear equations. 
It is provided, however, that the coefficient matrix of the simultaneous linear equations and 
the vector on the right part are more or less varied when the vibration is made with a 
sinusoidal wave, when the vibration is made with a random waveform, and when the input 
acceleration is made with a pulse waveform. The procedures will be shown below. 

(Mathematical 75) 



few 




Cl3 




Ct6 


Ct6 


<*,7 


Ct8 


Ct9' 










o& 


C Z3 








c Z7 


C 2,8 














C 3L3 


C&4 




Cae 




^33 










Ct2 




cm 


C 4JS 


c 46 


C4.7 






Syjg 










c w 








C*7 












c w 












<*7 














Gil 




C7.4 




C7.6 


C7.7 




<^.9 


S ajc 












<*4 






c ft7 


c ae 


C 8* 








1^1 




c a3 


<*4 


<*6 


^as 


C*7 


<?9jB 


<3Bj9; 









When the vibration is made with a sinusoidal waveform, a random waveform, or a 
pulse waveform, the coefficient C„ (I = 1 — 9, j =1 - 9 ), d; (i = 1 — 9) of the 
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aforementioned simultaneous linear equations is decided by the following procedures. 
1) The case of using sinusoidal waveform for input acceleration 
When the input acceleration is made with a sinusoidal waveform, theoretically the 

relation between the input acceleration vector (aix, a* y , aiz)exp(ja>t) and the output 

acceleration vector (ao X , ao y , ao Z )exp(ja>t) is as shown below. 
(Mathematical 76) 



Procedure 1 

Theoretically, an input acceleration is imparted first in the direction of vector (aixi, 
aiyi, aizi) at a given frequency fi(coi=27rfi). The input acceleration vector signal is 
represented by (aba, ai yJ , aizi)exp(jcoit). The X axis output signal, the Y axis output signal, 
and Z axis output signal of the accelerometer at this time are respectively represented by 
aoxi(a>i)exp(ja>it), aoyi(a>i)exp(ja>it) and ao 2 i(a>i)exp(j(Oit). 

Since the actual vibration is made with a burst signal and consequently the signal 
(ajxi, aiyi, aizOexpGoot) is also a burst signal, the three components ought to be considered 
as resulting from the extraction of the fi component through FFT. In other words, when 
the actual vibration signal is expressed on the time base by the function of the vector (aixi, 
aiyi, ai2i)Bi(t), the part Bi(t) means the burst signal. It is assumed that the complex 
number Z l& i results from extracting the fi component from the burst part through FFT. 

Since the output signals (B ox i(t), B oy i(t), B oz i(t)) of the accelerometer in this case 
are also burst signals, the fi component is extracted likewise through FFT. The fi 
component which is obtained at this time is a complex number, the division of this 
component by Zisi results in (aoxi(oi)i), aoyi(coi), aozsi(coi)). As a result, the following 
formula is established for the frequency component 1 . 



' a at exp(j(Dt) y 



Syvc S»y S y§z a^expijat) 
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(Mathematical 77) 
Syje &y.y 



V > 




'tad**? 















When this formula is unified as equations concerning (S w , Sx, y , S^z, S y , x , S y , y , S y>2 , 
Sz^, S^y, Sz, 2 ), the following formula is established. 
(Mathematical 78) 

C^=0 C^2=0 0^3=0 0^=0^ C2fi=a^ C^a^ 0^=0 ^=0 0^=0 

c*i=o c w =o c^a=o c w =o ^=0 c ze =o ^7=^ 0^=^ 0^9=^1 

rf 1=^al(«h) <*2=«<Brf(«l) *3=*«1<<»l) 



Procedure 2 

Theoretically, a vibration is imparted first in the direction of vector (aix2, ai y2 , aiz2) at 
a given frequency fi(coi=27tfi). The input acceleration vector signal is represented by (a^, 
ai y 2, aiz2)exp(j(0it). The X axis output signal, the Y axis output signal, and Z axis output 
signal of the accelerometer at this time are respectively represented by aox2(coi)exp(jcoit), 
ao y i(a>i)exp(ja>it) and ao Z 2(a>i)exp(jcoit). 

Since the actual vibration is made with a burst signal and consequently the signal 
(aba, ai y 2 3 aiz2)exp(ja>t) is also a burst signal, the three components ought to be considered 
as resulting from the extraction of the fi component through FFT. In other words, when 
the actual vibration signal is expressed on the time base by the function of the vector (a^, 
ai y 2, aiz2)B 2 (t), the part B 2 (t) means the burst signal. It is assumed that the complex 
number results from extracting the fi component from the burst part through FFT. 

Since the output signals (B ox2 (t), B oy2 (t), B oz2 (t)) of the accelerometer in this case 
are also burst signals, the fi component is extracted likewise through FFT. The fi 
component which is obtained at this time is a complex number, the division of this 
component by Z is2 results in (aox2(coi), ao y2 (ooi), aoz2(a>i)). As a result, the following formula 
is established for the frequency component 1 . 
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(Mathematical 79) 





S x,y 




( > 








s y.y 


S y.z 






«cy2(®l) 




S z.y 











When this formula is unified as equations concerning (S*^, Sx, y , Sx^, S y ^, S y>y , S y>z , 
Sz^c, Sz, y , S^z) to derive the coefficients of Mathematical 75, the following formula is 
established. 

(Mathematical 80) 

<?Vl=%2 Ci2 = q*£ Ci3=^£i2 C M = 0 C <5 =0 Gi6 = 0 Q,7=0 ^=0 C W =0 

C &1 = 0 C^ 2 = 0 Ce3=0 q BL4 =o 0^=0 c^ 6 = o c^ 7 = ^ c w = ^ C W=^2 
*4=*ar2<»1> ^=V^ d 6 = a oz2(*l) 

Procedure 3 

Theoretically, a vibration is imparted first in the direction of vector (aix3, ai y3 , a^) at 
a given frequency fi(a>i=27rfi). The input acceleration vector signal is represented by (a^, 
ai y 3, aiz3)exp(jcoit). The X axis output signal, the Y axis output signal, and Z axis output 
signal of the accelerometer at this time are respectively represented by aox3(a>i)exp(ja)it), 
aoy3(coi)exp(jcDit), ao23(coi)expGa>it). 

Since the actual vibration is made with a burst signal and consequently the signals of 
vibration vectors (a^, ai y3 , ai23)exp(ja)t) are also burst signals, the three components ought 
to be considered as resulting from the extraction of the fi component through FFT. In 
other words, when the actual vibration signal is expressed on the time base by the function 
of the vector (a^, ai y3 , aiz3)B 3 (t), the part B 3 (t) means the burst signal. It is assumed that 
the complex number Z is3 results from extracting the fi component from the burst part 
through FFT. 

Since the output signals (B OX 3(t), B oy3 (t), B oz3 (t)) of the accelerometer in this case 
are also burst signals, the fi component is extracted likewise through FFT. The fi 
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component which is obtained at this time is a complex number, the division of this 
component by Zui results in (aoxsfal), ao y3 (coi), aozsC^i). As a result, the following formula 
is established for the frequency component 1 . 
(Mathematical 81) 



When this formula is unified as equations concerning (Sx^, Sx, y , Sx, z , S y>x , S y , y , S y , z , 
Sz^c, Sz, y , Sz^), the following formula is established. 
(Mathematical 82) 

qM=o c w =o c w =o c M = ^3 c w =^a <^6 = ^3 ^7 = ° ^aa=° Qu> = ° 

Cfc,1=0 09,2=0 0^3 = 0 C M = 0 C M =0 qajs = 0 Cqj = <fe3 0^ = ^ 0^9 = ^ 
<h = *«3<<»|) = ^3^) d 9 = <fcr3<<»l> 

Procedure 4 

Consequently, the coefficient matrix of the aforementioned simultaneous linear 
equations has been decided and the vectors on the right part have been decided as well. 
The aforementioned simultaneous linear equations are solved without fail because the 
vectors in the three directions are selected in linear independence. Subsequently, the 
matrix sensitivity in a different frequency can be found by setting the vibration frequency at 
a different value and performing the same calculation. Specifically, the frequency is 
scanned in a necessary frequency band. 

The cases of using four-axis through six-axis accelerometer will be described by the 
following generalization for the sake of brevity of explanation. 

[Generalization] 

The degree of freedom of a vibration source as a motion generating machine is 
denoted by N and the number of axes of an inertia sensor by M. At this time, the 
following condition must be satisfied. 



(s 



'yjc S^ z 0^3 = acy3fc>l) 



S. 
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N < 6, M < 6 (Both do not become 1 simultaneously) 

The coordinate axes are denoted by Xi, X 2 , X 3 , X4, X 5 , X* and will be always 
enumerated in the order mentioned. The parameter of the N degree of freedom of a 
motion generating machine is denoted by j. Then, the number of axes of an inertia sensor 
is denoted by M and the parameter thereof by i. In this case, the following formulas are 
satisfied. 

1 < j < N, 1 < i < M 

At this time, the dynamic sensitivity matrix of an inertia sensor at a given frequency 
(co = 27tf) is expressed by a matrix of M x N (M rows and N columns). Let N stand for 
the degree of freedom of a motion generating machine and M for the number of axes of an 
inertia sensor, and it will become possible to decide the matrix of the order of M x N as the 
sensitivity. 

(Mathematical 83) 

f &n(fl>) S\jf(a) y 

* * 

The simultaneous linear equations for deciding all the components of the sensitivity 
matrix are given by any simultaneous linear equations other than those which have MN 
lines and MN columns of matrix as a coefficient matrix and constant vectors having MN 
components as a constant term. 
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(Mathematical 84) 



r c„ 



^2,2 



c 



V 



^ MN,MN ) 



4 1 



4* 



The procedure for deciding the coefficient matrix and the constant term of the 
simultaneous liner equations (Mathematical 84) more or less varies when the input 
acceleration is made with a sinusoidal wave burst waveform, when the input acceleration is 
made by a random vibration, or when the input acceleration is made with a pulse waveform. 
The procedures which are contemplated here will be described below in the order of steps. 

[The case of using a sinusoidal wave burst signal for the input acceleration] 

Since the vector space which can be formed by a vector space generating machine 
has N dimensions, N linearly independent vectors can be selected. So, the vibration 
vector, the burst signal pertinent then, and the output signal are applicable as shown in the 
following table. 

(Table 14) 



N oscillation 
vectors 


Oscillation vector 


6) j component 
resulting from 
spectral analysis 
of burst signal 


Output signal of inertia 
sensor 


Number of angular 
frequency resulting from 
spectral analysis of 
output signal: a>! 
component 


First vector 




Z\> hi, l 




Zb.bs.xl. 1.1 * 

*"» ^j.bs.xM.1.1 


Second vector 




bs.2 


■ "W 2 (t)) 


Zb,bi. xl.Zl> 

•",Z b .b».xM.2.1 


Third vector 


-Ol.il a) Bs© 


bo 




^».bs.xi.ii» 

Zb.bi.aM. 3. 1 












N'th vector 


( a txl.H» a LiZH> 

-airiiH)BH(t) 




(WHO), ao.xzn(t), 

•••^.,M.N(t)) 


Zb.bs. al.N.l> 

" '.ZbJa.rM.N.l 
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The coefficients of the simultaneous linear equations of (Mathematical 84) are 
selected from this table as follows. 

Procedure 1 (which corresponds to the first vibration vector) 
(Mathematical 85) 

■ 

*\ 53 Z atejrtll/ Z Uu.1 = Z aiu2.ll/^fe1 - - - <Sf = Z «*u*f.tl/ Z tM 

Procedure 2 (which corresponds to the second vibration vector) 
(Mathematical 86) 

m 

m 

Procedure 3 (which corresponds to the third vibration vector) 
(Mathematical 87) 

4a*+1 = z afcu1Al/^.6^3 *2ar+2 = ^aftwZai/ 2 ^ ^ = z *fa.*M3,l/ Z a*3 
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Procedure N (which corresponds to the Nth vibration vector) 
(Mathematical 88) 



^-1)1/42.1 = C (JT-t)Jf+2£ = ^-11^+2^ ~ 0 

■ 

C A0kf.1 — ' 83 C XMiM-W m 0 

<^1)*f+1 = Z aJbs*\NAl Z lJxJf &1 = Z oJbt^Zw/ Z iMM rf AH/ = ^fisjAdjn/^MN 

Procedure N + 1 

By the preceding procedures, all the components of the coefficient matrix and the 
constant terms of the equations (Mathematical 84) have been decided and the fact that the 
determinant of the coefficient matrix is not zeroed when the vibration vector is linearly 
independent can be proved. The solution aimed at, therefore, is attained without fail. 
Thus, all the coefficients of the matrix sensitivity are uniquely found. What is to be done 
resides in having the frequency of vibration scanned as occasion demands. 

[The case of using random signal for the input acceleration] 

Since the vector space which a motion vector generating machine can generate has 
N dimensions, N linearly independent vectors can be selected. So, the vibration vector, 
the results of the spectral analysis of the random signal pertinent thereto, and the output 
signal are applicable as shown in the following table. In short, due attention should be 
paid to the fact that the expression "random vibration" does not necessarily imply that the 
vibration is random also in direction. It requires attention that the frequency component is 
random and the N linearly independent directions are fixed. At this time, the following 
table is obtained. 
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(Table 15) 



N oscillation 
vectors 


Oscillation vector 


CO] component 

ramltina Frnm 

resuiiuig num 
spectral analysis 

of random 
oscillation signal 


Output signal of inertia 
sensor 


Number of angular 
frequency resulting from 
spectral analysis of 
output signal: CO] 
component 


First vector 






(a^xi.iCt), ^^(t), 

"Witf)) 


Zb,bi.xl.l.l» 

Zb.bs.xM. LI 


Second vector 


0 3 i.*1.2»< 3 i.x2.2> 

•••ai.xN.2)B2(t) 


bs. 2 


(So. xl^(t), a^C), 


^.bs.xM.2.1 


Third vector 


• aixN.3) BgCt) 




(ao.n. 3 (t), ^^(t), 


^bs. sl.3,l> 












]STth vector 


" a ixRK ) B^t) 




(ao,xi,N(t), kAMv). 


^»,bs.xl.N,l> 



Since the coefficient matrix and the constant term vector of the simultaneous liner 
equations (Mathematical 84) are decided by carrying out a procedure similar to that 
adopted in the case of a sinusoidal wave, all the components of the matrix sensitivity 
regarding the a>i component can be decided. When the resolution found by the spectral 
analysis of a random signal is Ago, all the components of the matrix sensitivity at necessary 
angular frequencies can be found by selecting the angular frequency of the formula of 
(coi ± Aa> x an integer) and following the same procedure. This is because the vibration 
vector is linearly independent and the coefficient determinant of the simultaneous linear 
equations (Mathematical 84) is not zero at any angular frequencies. 

[The case of using a pulse waveform signal for the input acceleration] 
Since the vector space which can be formed by a motion generating machine has N 
dimensions, the N linearly independent vectors can be selected. Thus, the spectral analysis 
of a pulse signal which exists at the vibration spectrum has shown that the output signal, for 
example, can be written as shown in the following table. At this time, the following table 
is derived. 



(Table 16) 
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N oscillation 
vectors 


Oscillation vector 


<Oj component 
resulting from 

spectral analysis 
of pulse 

oscillation signal 


Output signal of inertia 
sensor 


Number of angular 
frequency resulting from 
spectral analysis of 
output signal: CO ! 


First vector 




^ hs.l 


(a-.xi.i(t).^ ia . 1 G), 


component 


Second vector 


( fl ixl.2><2i.xZ2> 


2 


(a D .x 1 .2(t),a D . ia2 (t), 


"y Z^bs. ^ j ! 

^b, bs. xl.2,l> 


Third vector 






,fttt 


'"> ^ta.hi.iHyr.^i 










, ""> ■ 


>Tth vector 








» ^b.t*. xM, N, 1 



Subsequently, by following the same procedure as in the case of a sinusoidal wave 
« .a made possible to decide the coefficient matrix of the simultaneous linear equations 
(Mathematical 84) and the constant term vector and consequently decide all the 
components of the matrix sensitivity regarding the angular frequency <o, When the 
resolution found by the spectral analysis of a pulse signal is Ao>, all the components of the 
matnx sensitivity a, necessary angular frequencies can be found by selecting the angular 
frequency of «he formula of (o,±A«, x an integer) and following the same procedure 
Thta ,s because the vibration vector is linearly independent and the coefficient matrix of the 
stmultaneous linear equations (Mathematical 84) is no. zero a, any angular frequencies 

According ,o the first aspect of tins invention, by causing a plurality 0 f hser 
mterferometers disposed aiong as many directions to seize multidimensional motions 
entenng an ntertia sensor as input signals and comparing these input signals with output 
stgnals from the inertia sensor, the sensitivity of the inertia sensor cat, be defined in the 
rnatrtx form. The result is tha. tins inertia sensor has been calibrated. While the 
conventional captation technique has enabled the accelerometer to be calibrated solely 
wth the amplitude of the acceleration, the technique of this invention enables the 
accelerometer to be calibrated by the acceleration as a vector and enables the accelerometer 
to measure the acceleration. 
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Then, according to the second aspect of this invention, the sensitivity is measured in 
the matrix form under an epochal input acceleration waveform and the calibration can be 
carried out in this framework. 

According to the third aspect of this invention, the sensitivity is measured in the 
matrix form under an input acceleration waveform of the nature of the pulse function and 
the calibration can be carried out in this framework. 

According to the fourth aspect of this invention, the sensitivity relying on the 
frequency is measured in the matrix form under the conditions of the input acceleration of 
the nature of the pulse function and the calibration is carried out in this framework. 

According to the fifth aspect of this invention, the sensitivity is measured in the 
matrix form under the random input acceleration waveform and the calibration can be 
carried out in this framework. 

According to the sixth aspect of this invention, the sensitivity relying on frequency 
is measured in the matrix form under the random input acceleration condition and the 
calibration can be carried out in this framework. 

According to the seventh aspect of this invention, the measurement can be fulfilled 
by the use of a laser interferometer while the motion generating machine is continuously 
monitored. 

According to the eighth aspect of this invention, the aforementioned monitoring can 
be carried out with high accuracy. 

According to the ninth aspect of this invention, the error can be evaluated. 
According to the 10 th aspect of this invention, when the laser interferometer is not 
easily used and some other alternative device of high precision is adopted, the regression as 
to the precision of measurement can be retained to a slight extent. 

According to the 1 1 th aspect of this invention, the influence of the gravitational 
acceleration can be excluded from the measured value. 

According to the 12* aspect of this invention, it is made possible to define and 
observe the cross or transverse sensitivity easily. While the piezoelectric type one-axis 
accelerometer heretofore widely used as .the standard accelerometer has failed to afford 
cross or transverse sensitivity as complex numbers, this invention affords the cross or 
transverse sensitivity as complex numbers. 
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According to the 13 th aspect of this invention, the reliability of the measurement 
with the accelerometer which is used in an applied field in which the inability to make 
three-dimensional measurement as in the seismometer used for measuring an earthquake 
comprising pitching and rolling proves meaningless can be exalted. 

According to the 14 th aspect of this invention, the disturbance to the vibro-isolating 
common bed is detected more accurately and the effect of vibration insulation is enhanced. 
The application of the invention to the vibro-isolating device, for example, results in 
contributing to the promotion of miniaturization of the semiconductor. The realization of 
the evaluation of performance of the multi-degree freedom inertia sensor (particularly made 
of semiconductor) simultaneously possessing translational motion acceleration detecting 
function, angular velocity detecting function, and angular acceleration detecting function 
results in improving the precision of manufacture of a semiconductor device and exalting 
the performance of the product. 

According to the 15 th aspect of this invention, the technique of utilizing the 
multi-axis accelerometer which is used for controlling the suspension of an automobile or 
controlling the automobile passenger protecting airbag can be exalted. 

According to the 16 th aspect of this invention, the precision of the control of the 
movement of a robot is exalted. In the mechanism such as a man-like robot which 
generates an unrestricted motion in conformity with the situation inside a three-dimensional 
space while it is continuously subject to a feedback and control, the use of the inertia sensor 
calibrated with the vector acceleration results in enabling the motion to be generated in 
greater fineness than ever. 

According to the 17 th aspect of this invention, the application such as the 
accelerometer installed in the interior of an automobile for ensuring safety of the passenger 
in the case of a car collision in which the inability to measure the acceleration as the vector 
proves meaningless enjoys an increased utilization and ensures a heightened safety of an 
automobile. 

According to the 18 th aspect of this invention, even an acceleration measuring 
device, an angular velocity sp e ed measuring device, or an angular acceleration measuring 
device which is furnished with an output only in one axis is enabled to discern the 
sensitivity in other axis. 
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According to the 19 th aspect of this invention, the inertia sensor is enabled to effect 
the measurement of dynamic matrix sensitivity even when the calibration is required to be 
effected with a greater degree of freedom than the degree of freedom of the motion 
possessed by the motion generating machine or when the number of axes which can be 
handled is smaller than the number of axes required to be used in the calibration. 

According to other aspects of this invention, the sculling error, i.e. an error due to 
the rotation of an accelerometer, which has been heretofore measured with difficulty, can 
be measured. The term "sculling error" as used herein means an error ascribable to an 
apparent acceleration (rectifying acceleration) which arises when the linear vibration 
acceleration generated along a given axis and the angular acceleration having the same 
frequency as the linear vibration acceleration and generated around an axis perpendicularly 
intersecting the given axis occur simultaneously. This acceleration has a magnitude which 
depends on the amplitude and the phase of the relevant inputs and appears in an axis 
perpendicular to the two axes. 

Industrial Applicability: 

This invention is configured as described above. By using the device of this 
invention for measuring the dynamic matrix sensitivity of an inertia sensor and the method 
of this invention for operating the device, therefore, it is made possible to acquire exactly 
the multidimensional sensitivity property of an inertia sensor, judge whether or not a given 
inertia sensor fits the purpose of use, and provide an accurate data for designing an inertia 
sensor which fits the purpose of use. It is further made possible to design and construct a 
still more complicated complex inertia sensor by making use of the inertia sensor which has 
secured the sensitivity property accurately. 



